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COMPUTING PROGRAM FOR AXIAL DISTRIBUTION OF AERODYNAMIC 
NORMAL- FORCE CHARACTERISTICS FOR AXISYMMETRIC 
MULTISTAGE LAUNCH VEHICLES 

By Ragan B. Madden 
Langley Research Center 

SUMMARY 

This report describes a digital computer program which calculates the axial distri- 
bution of aerodynamic normal-force characteristics for axisymmetric, multistage launch 
vehicles in the linear angle-of-attack range. In calculating these characteristics for a 
particular configuration the program utilizes available experimental data, the majority of 
which were obtained from NASA Technical Note D-3283 and Lockheed Missiles & Space 
Company report LMSC/HREC A712618. These data cover a wide range of Mach numbers 
and geometric parameters and can be readily utilized to obtain similar characteristics 
for most launch- vehicle configurations. 

INTRODUCTION 


The increased demand upon the aeroelastician to obtain the distributed aerodynamic 
normal-force loads for launch- vehicle configurations necessitated the development of a 
digital computer program to calculate these force coefficients. 

A large amount of experimental data for the components comprising most launch 
vehicles were obtained from references 1 and 2. These data were reduced to a standard 
form and employed in the program as fixed input data. In obtaining the characteristics 
for a particular vehicle, the program considers the individual components which comprise 
the vehicle and performs linear interpolations or extrapolations on the proper fixed input 
data values. 

The methods of reducing the data to a standard form and the limitations on these 
methods are described in reference 1, which should be used in conjunction with this 
program. 

This report consists of two main parts. The first part is the users' section which 
is to be referred to when one wants to obtain the characteristics for a particular vehicle. 
This part describes the input and the output, and two examples are given in order to illus- 
trate the use of the program. The second part is the programer's section which describes 



the program logic and the method of inputting the fixed data. If one is just interested in 
using the program to obtain the characteristics for a particular vehicle, then this part of 
the report need not be referred to. Also included are the program listing and the fixed 
input data listing as appendixes A and B, respectively. 

SYMBOLS 


When components are mentioned in series, they are to be considered in sequence 
from nose to aft end. 


A 

A1 


experimentally determined coefficient to determine X values for a sphere 


denotes values of interpolated X values at (fa)r anc * a *- eac ^ 

value in figure 10; denotes values of interpolated X values at and 

q)k2 


VU.J.UV in ugui 

at each values in figure 11 


A2 denotes interpolated X values at (9f ) r and at each ( value in fig 

/Di\ ^ D °^ 2 /x\ 

ure 10; denotes interpolated X values at I— and at each [— ) 

\pJr \ D o/k2 


value in figure 11 


A3 


denotes interpolated X values at (0 n )r and at each value i 

figure 10 ' oJ\& 


A4 


A44 


denotes final interpolated X values at M r and at each f-2-j value 

\P°/k2 

denotes interpolated X values at (0 n ) r and at each value in 

figure 9 k2 


C N normal-force coefficient, F^qS 


C N, 


a 


normal-force-coefficient slope, , per radian 

da 


D 


local diameter, meters 
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f 


Do 

D 1 

d 2 

d 3 

d 4 

DX4 

DX5 

DX7 

( D l/°o)i 

DX 

Fn 

f a 

f n 

L 


reference diameter, meters 

reference diameter of first frustum, meters 

reference diameter of second frustum in multistage vehicle, meters 

reference diameter of third frustum in multistage vehicle, meters 

reference diameter of fourth frustum in multistage vehicle, meters 

diameter of second cylinder of ogive-cylinder-cylinder combination, meters 

diameter of second cylinder of a cone-cylinder-cylinder combination, meters 

diameter of second cylinder of a frustum-cylinder-cylinder combination, 
meters 

reference diameter of ith frustum divided by reference diameter of cylinder 
following ith frustum 

product of local diameter and generalized loading function, referred to in this 
report as load characteristics, meters per radian 

normal force, newtons 

fineness ratio of cylinder preceding frustum 

fineness ratio of tangent ogive 

length of component under consideration, meters 
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M mach number 

l length of sharp ogive preceding cylinder or length of cylinder preceding 

frustum, meters 

q dynamic pressure, newtons per meter2 

R planview radius of ogive, meters (see fig. 5(b)) 

r one-half reference diameter of ogive, meters 

S reference area, meters2 

x local body axial coordinate of component, p - p Q , meters 

XI junction point of cylinder and following frustum, meters 

X2 junction point of frustum and following cylinder, meters 

XR1C junction point of cone-cylinder or ogive- cylinder, meters 

XR4 p coordinate at which second cylinder of ogive-cylinder-cylinder combination 

begins, meters 

XR5 p coordinate at which second cylinder of cone-cylinder-cylinder combination 

begins, meters 

XR7 p coordinate at which second cylinder of frustum -cylinder -cylinder combi- 

nation begins, meters 


4 



dCjr 

S - product of reference area and distributed normal-force-coefficient slope, 

dx 

meters per radian 


a 


angle of attack, radians 


e b 


boattail angle, degrees 


0 f,i 


ith frustum angle where, in multistage vehicle, ith frustum is defined to be 
closer to nose than i + 1st frustum 


cone semivertex angle, degrees 


1 dCN 0! 

generalized loading function, — S , per radian 

D dx 


axial coordinate of launch vehicle (or configuration), meters 


fJ-o 


value of p at origin of component under consideration such that locally 
x = p - fi 0 , meters 


Subscripts : 

i,j,k denote fixed input variables in figures 9, 10, and 11 

k2 ranges over number of case input x stations for a given component 


n nose 


r case input variable in figures 9, 10, and 11 
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T 


tangency point 


1,2 denote the first and second interpolations, respectively, at a given case input 

variable in figures 9, 10, and 11 

APPLICATIONS 

A comparison of load characteristics from experimental and empirical data for a 
typical launch vehicle is shown in figure 1. Empirical load characteristics for typical 
launch vehicles are shown in figure 2. 



Figure 1.- Comparison of load characteristics from experimental and empirical data for typical launch-vehicle configuration. 
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USERS' SECTION 


CONFIGURATIONS 


From references 1, 2, and 3, X values for eight basic components were obtained. 
These values are read into the program as fixed input data each time the program is used. 
The description of the components, the component code numbers, and the parameter and 
Mach number ranges covered by these X values are shown in table I. The solid-line 
components in table I are the components of interest, and the dashed-line components are 
the components which precede the components of interest. 


TABLE I.- COMPONENTS FOR WHICH LOADING FUNCTIONS (X VALUES) ARE AVAILABLE 
AND PARAMETER AND MACH NUMBER RANGES COVERED 


I 

Cone - sharp or blunted 


Tangent ogive - sharp or blunted 


3 

Sphere 


Cylinder following sharp or 
blunted tangent ogive 


Cylinder following sharp or 
blunted cone 


Frustum following sharp or 
blunted cone - cylinder 


7 

- r -r ^ 

I T 1 

t 

Cylinder following cone- 

< l D f ! 

Do 

cylinder-frustum 


_L 


Boattail following cylinder 
at aft end 


Diagram of component 


~T 

I 

D 0 

A 




k l h 


STT 

\ D 0 


0f,l 


"tf" 


- I - 


T" 

Do 


^b 


I D 0 


parameter range 


0° £ fl £ 50° 


0.50 £(f n = L/D 0 ) = 6.0 


3s(f n = l/D o )?7.0 


0° ^ 0. i 40° 


15° £ e n S 30° 

o° 5 e if 1 § 2o° 


0° £ 0 f ! S 20° 

0 5 (D^DoJi S 1 


Mach number range 


0.70 to 15.00 


0.80 to 15.00 


0.25 to 15.00 


3.00 to 15.00 


0.70 to 15.00 


0.80 to 15.00 


0.80 to 15.00 


1.50 to 15.00 
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Figure 3 shows the geometric parameters which are used in defining the eight basic 
components, and figure 4 shows five variations of the eight basic components for which 
the program may be used. The solid and dashed lines in figures 3 and 4 have the same 
significance as in table I. 



Sharp or blunted 
cone 



Sharp or blunted 
ogive 



Spherical segment Cylinder following sharp or blunted 

tangent ogive 




Cylinder following sharp 
or blunted cone 



Frustum following cone- 
cylinder 


Figure 3.- Parameters used to define geometrical components. 
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I 


I 

I 




Boattail following cylinder at aft end 



(D../D = D-VD 

v V o'l Vo 


Cylinder following cone-eylin&er-frustum 
Figure 3.- Concluded. 




Station XR5 


Variation I: Cylinder -cylinder following cone 



Station XR4 


Variation II: Cylinder -cylinder following ogive 



Station XR7 


Variation III: Cylinder-cylinder following frustum 



Variation IV: Frustum following ogive -cylinder 

(Restricted to f a = 7/D 0 - 4.0) 



Variation V: Multifrustum-cylinder combinations 

Figure 4.- Five variations of components 4, 5, 6, and 7 and associated parameters. 
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CASE INPUT 


This section gives the method of inputting data for a particular vehicle to the pro- 
gram. Table n gives the FORTRAN names of the input variables and the definitions of 
each variable. 


TABLE H.- FORTRAN NAMES AND DEFINITIONS FOR CASE INPUT DATA 


FORTRAN name 

Definition 

COMPNO(I) 

Component code numbers, beginning with component at nose, I = 1 to NCMPS; 
see table I and restrictions (5) and (6) 

DDRR(J) 

(Di/Do)^ J = 1 to NFCYCM; the reference diameter of Jth frustum divided by 
reference diameter of cylinder following Jth frustum; see figures 3 and 4 

DO(I) 

Reference diameter Do of each component, beginning with component at nose, 
I = 1 to NCMPS; see figure 3 

DX4 

DX4, diameter of second cylinder of ogive -cylinder -cylinder combination; see 
figure 4 

DX5 

DX5, diameter of second cylinder of cone- cylinder- cylinder combination; see 
figure 4 

DX7 

Diameter of second cylinder of frustum- cylinder- cylinder combination; see 
figure 4 

FAR 

f a , fineness ratio of cylinder preceding frustum; see figure 3 and restriction (8) 

FNR 

f n , fineness ratio of sharp tangent ogive; see figure 3 and restriction (13) 

MACHN(K) 

Mach numbers at which DA values are desired, K = 1 to NMACHN; see 
restriction (1) 

NCI 

Input as 1 if variation I exists and input as 0 otherwise; see figure 4 and restric- 
tion (6) 

NC2 

Input as 1 if variation II exists and input as 0 otherwise; see figure 4 and restric- 
tion (6) 

NC3 

Input as 1 if variation m exists and input as 0 otherwise; see figure 4 and restric- 
tion (6) 

NCMPS 

The number of components which comprise the vehicle; see restrictions (5) and (6) 

NDX 

Input as 0 if stations are to be input from nose to aft end in ascending order and 
input as 1 if stations are to be input from nose to aft end in descending order 

NFCYCM 

The number of frustum- cylinder combinations; see restriction (7) 

NMACHN 

The number of Mach numbers; see restriction (1) 

NOXS 

The number of stations; see restriction (2) 

NPUNCH 1 

Input as 0 if no punch card output is desired and input as 1 if punch card output 
is desired 

NTHFRR 

The number of frustum angles 0f; see restriction (7) 

NXS(I) j 

The number of x stations of component I, beginning with component at nose, 
I = 1 to NCMPS; see restriction (4) 

THBR 

Q b, boattail angle in degrees; see figure 3 

THFRR(N) 

Of n , Nth frustum angle in degrees, beginning with frustum closest to nose, 
N = 1 to NTHFRR; see figures 3 and 4 and restriction (7) 

THNR 

6 n , cone semivertex angle in degrees; see figure 3 

XR1C 

XR1C, junction station of cone- cylinder or ogive -cylinder; see figure 3 

XR4 

XR4, fj. coordinate at which second cylinder of an ogive-cylinder-cylinder com- 
bination begins; see figure 4 

XR5 

XR5, /lx coordinate at which second cylinder of a cone -cylinder- cylinder begins; 
see figure 4 

XR7 

XR7, [i coordinate at which second cylinder of a frustum -cylinder -cylinder 
begins; see figure 4 

XS(M) 

p. coordinates at which DA values are desired, beginning with station at nose, 
M = 1 to NOXS; see restrictions (2) and (3) 



Restrictions on Case Input 


The following restrictions apply to the case input: 

A maximum of 15 Mach numbers may be input. 

A maximum of 350 coordinates overall may be input, and they must be input from the 
nose to the aft end in either ascending or descending order. These stations are the 
XS(I) values defined in table n. 

At junction points of components — excluding variations I, n, and in - the associated 
stations must be input twice in the XS(I) set of case input. At junction points of 
variations I, H, and HI - that is, at stations XR4, XR5, and XR7 - the stations are 
input once on the fourth case input data card and once in the XS(I) set of case input. 

A maximum of 100 stations for an individual component may be input. 

The component code numbers COMPNO(I) refer only to the eight basic components 
shown in table I. It should be understood that a blunted cone is composed of compo- 
nents 3 and 1 and a blunted ogive is composed of components 3 and 2. The cylinder- 
cylinder combination for variations I, II, and DI is treated as a single component - 
that is, for variation I, the cylinder -cylinder combination is component 5; for varia- 
tion II, component 4; and for variation m, component 7. The presence of variations I, 
n, or m in a vehicle is denoted by the value of the case input variables NCI, NC2, 
and NC3, respectively, on the second case data card. Variation IV is just a particu- 
lar arrangement of some of the eight basic components. The presence of variation V 
in a vehicle is denoted by the value of the case input variable NFCYCM on the sec- 
ond case data card. The values of COMPNO(I) must always identify the components 
taken in sequence from the nose to the aft end, and the value of NCMPS must be 
less than 15. 

A vehicle is not allowed to be composed of two or more variations of the same type. 

As an example, a vehicle is allowed to have one variation I and one variation n but is 
not allowed to have two variations I. 

The number of frustum -cylinder combinations may not exceed five. If there are no 
such combinations or only one such combination, the value of the case input variable 
NFCYCM must be 1. Also, the number of frustums may not exceed five. If these 
are no frustums or only one frustum, the value of the case input variable NTHFRR 
must be 1. 

If more than one frustum-cylinder combination is present, the value of the case input 
variable FAR is the fineness ratio of the cylinder preceding the frustum closest to 
the nose. If the first frustum follows variations I, II, or HI, then the input value of 
FAR is given by 



_ Length of l a rger diameter cylinder + Length of smaller diameter cylinder 
Diameter of larger diameter cylinder 

(9) For variations I, n, and m, the vehicle cannot have a third cylinder immediately fol- 
lowing the smaller diameter cylinder. This cylinder can be followed only by a boat- 
tail, frustum, or frustum -cylinder combination for variations I and II and only by a 
boattail for variation in. 

(10) For variation V a variation in is allowed only after the last frustum-cylinder 
combination. 

(11) For variation IV, the value of the case input variable THNR must be 0.0. 

(12) The units of the input station coordinates and diameters are optional but should be 
self-consistent. The units of the input angles must be degrees. 

(13) The value of the case input variable FNR must be the fineness ratio of the sharp 
tangent ogive. If the ogive is blunted, then the length of the sharp ogive can be 
determined by 


L = 



and the value of FNR is L/D 0 . 

Considerations Which Should Be Given to Case Input 

It is desirable to input the stations for each component at approximately every 0.2D o 
for 0 ^ x/Do = 1.0, at every 0.25D o for 1.0 < x/D 0 = 3.0, and at every 0.5D o for 
3.0 < x/D 0 . Although this convention is not required, following it will generally yield an 
adequate description of the DA output. 

Data for more than one vehicle may be input at a time simply by placing the data cards 
for the second vehicle immediately after the data cards for the first vehicle. 

To obtain the DA distributions for a given vehicle over a range of Mach numbers, 
the outlined procedure should be followed. 

(1) The first step is to see if the program can be used for the vehicle. Table I and 
figure 4 show the components and variations for which the program can be used; hence the 
user must determine whether the components which comprise the vehicle are included in 
this table or figure. If they are, then the program can be used. If, however, a vehicle 
has a component which is not included in table I or figure 4, then the program may be used 
to obtain DA distributions up to this component but not aft of it. A common example of 
this type vehicle is one in which a boattail occurs in the middle of the vehicle - compo- 
nent 8 is a boattail at the aft end of a vehicle. In this case the program can obtain DA 
distributions up to the boattail but not aft of it. 
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(2) The section "Restrictions on Case Input" should be reviewed. 

(3) The case data must be properly coded and punched. This process is described in 
detail in the subsequent section "Case Input Data Cards Description." 

(4) The case input data cards must be placed after the last fixed input data card. 

(5) If DA values for more than one vehicle are desired, then the case input data for 
the second vehicle are placed immediately after the case data for the first vehicle. 

(6) The fixed input data must be input into the program each time the program is run. 
If, however, DA values for more then one vehicle are desired, the fixed input data are 
input only once. 


Case Input Data Cards Description 

The case input data are input on punched cards and, for each vehicle, the case data 
cards consist of 8 single cards followed by three sets of cards where each set contains 
from 1 to 35 cards. The case data variables are defined in table n. The values of these 
variables are either integer numbers (numbers punched without a decimal) or decimal 
numbers (numbers punched with a decimal) . Each integer number must be punched such 
that the last digit ends in a column which is a multiple of 5. The FORTRAN format used 
for the integer numbers is 1415, and no more than 14 numbers may be punched on a card. 

The decimal numbers must be punched according to the following restrictions. 

(1) If the case input decimal number is a whole number, 27 for example, then this num- 
ber may be punched either as 27. in which case the decimal point must end in a column 
which is a multiple of 7 or as 27.0 in which case the 0 must end in a column which is a 
multiple of 7. 

(2) If the number is positive, then the sign does not have to be punched. If the number 
is negative, the sign must be punched. 

(3) The number of values to the right of the decimal point may vary from 0 as in 27. 
to 6 as in .000001 but the total number of columns punched including the sign (if negative) 
and decimal point for a given number may not exceed 7. For example, the number 
-1580.12 cannot be punched as it is. It would have to be rounded off to -1580.1 in order to 
be punched. 

(4) The FORTRAN format used for the decimal number is 10F7.3, and no more than 
10 numbers may be punched on a card. 

The examples of the case data given in the following description are applicable for 
example configuration I which is shown at the end of the Users* Section in figures 5 and 6 
and table m. 
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First data card .- The first data card consists of a short heading which gives the name 
of the vehicle and any other information the user desires. The characters punched on 
this card are alphameric, and the first character is punched in column 2 and the last 
character cannot be punched beyond column 67. Spaces may be skipped between 
characters . 

Example : 


example: configukati on i 

Second data card .- The second data card consists of the values of the case input vari- 
ables NCMPS, NMACHN, NOXS, NDX, NCI, NC2, NC3, NFCYCM, NTHFRR, 
and NPUNCH. These values must be integers, and each value must end in a column 
which is a multiple of 5; that is, the value of NCMPS must end in column 5, the value of 
NMACHN must end in column 10, and so forth. 

Example : 

5 2 44 1 001 1 1 0 

Third data card .- The third data card consists of the values of the case input varia- 
bles THNR, THBR, FAR, and FNR. These values must be decimal numbers, and each 
value must end in a column which is a multiple of 7; that is, the value of THNR ends in 
column 7, the value of THBR ends in column 14, and so forth. Furthermore, there must 
be a value punched for each variable. If the vehicle does not have one of these variables, 
then the value of this variable must be punched as 0.0. 

Example : 


20.0 0.0 1.0 0.0 

Fourth data card .- The fourth data card consists of the values of the case input varia- 
bles XR4, XR5, XR7, DX4, DX5, DX7, and XR1C. These values are decimal num- 
bers, and each value must end in a column which is a multiple of 7; that is, the value of 
XR4 must end in column 7, the value of XR5 must end in column 14, and so forth. A 
value must be punched for each variable. If one of these variables is not applicable for 
the vehicle, then its value must be punched as 0.0. 

Example : 


0.0 0.0 25.0 0.0 0.0 8.0 165.0 


17 


Fifth data card .- The fifth data card consists of the values of the case input variables 
THFRR(I). The number of values of this variable is equal to the value of the case input 
variable NTHFRR punched on the second data card. The values of THFRR(I) are deci- 
mal numbers, and each value must end in a column which is a multiple of 7; that is, the 
value of THFRR(l) must end in column 7, the value of THFRR(2) must end in 14, and 
so forth. 

Example : 


12.5 

Sixth data card .- The sixth data card consists of the values of the case input variables 
DDRR(I). The number of values of this variable is equal to the value of the case input 
variable NFCYCM punched on the second data card. The values of DDRR(I) are deci- 
mal numbers, and each value must end in a column which is a multiple of 7; that is, the 
value of DDRR(l) must end in column 7, the value of DDRR(2) must end in column 14, 
and so forth. 

Example : 


.759 


Seventh data card .- The seventh data card consists of the values of the case input 
variables COMPNO(I). The number of values of this variable is equal to the value 
of the case input variable NCMPS punched on the second data card. The values of 
COMPNO(I) are integers, and each value must end in a column which is a multiple of 5; 
that is, the value of COMPNO(l) must end in column 5, the value of COMPNO(2) must 
end in column 10, and so forth. 

Example : 


3 15 6 7 

Eighth data card .- The eighth data card consists of the values of the case input vari- 
able NXS(I). The number of values of this variable is equal to the value of the case input 
variable NCMPS punched on the second data card. The sum of these values must equal 
the value of the case input variable NOXS punched on the second data card. The values 
of NXS(I) are integers, and each value must end in a column which is a multiple of 5; 
that is, the value of NXS(l) must end in column 5, the value of NXS(2) must end in 
column 10, and so forth. 
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Example: 


3 5 8 7 21 

First set of data cards follo wi ng first eight single data cards .- This first set of data 
cards consists of the values of the case input variable DO(I). The number of values of 
this variable is equal to the value of the case input variable NCMPS punched on the 
second data card. The values of DO(I) are decimal numbers, and each value must end 
in a column which is a multiple of 7; that is, the value of DO(l) must end in column 7, 
the value of DO(2) must end in column 14, and so forth. If the value of NCMPS is 
larger than 10 then the value of DO(ll) must end in column 7 of the next card, and so 
forth. 


Example : 


6.0 6.735 6.735 6.735 11.513 

Second set of data cards following first eight single data cards .- This second set of 
data cards consists of the values of the case input variable MACHN(I). The number of 
values of this variable is equal to the value of the case input variable NMACHN punched 
on the second data card. The values of MACHN(I) are decimal numbers, and each value 
must end in a column which is a multiple of 7; that is, the value of MACHN(l) must end 
in column 7, the value of MACHN(2) must end in column 14, and so forth. If the value 
of NMACHN is larger than 10, then the value of MACHN(ll) must end in column 7 of 
the next card, and so forth. 

Example : 


1.0 2.0 


Last set of data cards.- This last set of data cards consists of the values of the case 
input variable XS(I) . The number of values of this variable is equal to the value of the 
case input variable NOXS punched on the second data card. The values of XS(I) are 
decimal numbers, and each value must end in a column which is a multiple of 7; that is, 
the value of XS(1) must end in column 7, the value of XS(2) must end in column 14, 
and so forth. If the value of NOXS is larger than 10 (which is almost always the case), 
then the value of XS(ll) must end in column 7 of the next card, and so forth. 
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Example : 


175.0 

174.0 

173.0 

173*0 

171.0 

169.0 

1 67.0 

1 65.0 

165.0 

164.0 

163.0 

162.0 

161.0 

160.0 

15 b . 01 56 . 265 156. 265 

1 55.0 

1 54 . 0 

153.0 

152.0 

151.0 

1 50.0 

150.0 

14 b . 0 

146.0 

144.0 

142.0 

140*0 

138.0 

1 36.0 

134.0 

1 30.0 

125.0 

120.0 

115.0 

. 110.0 

1 00.0 

75.0 

50*0 


25.0 15.0 5.0 - 5.0 

OUTPUT 

The output begins with a listing of the input parameters and code numbers, and, for 
each Mach number, the stations and load characteristics (DA values) follow. For ogive- 
cylinder, sphere-cone, and sphere-ogive combinations, the DA values are continuous at 
the junction points and only one value is printed at the junction points. For each of the 
other combinations, the DA values are discontinuous at the junction points and two values 
are printed at each junction point. 

For components 4 and 8 no fixed input data in the subsonic and transonic range are 
input to the program. For any case input Mach number which is less than or equal to 1.35, 
the message NO FIXED INPUT DATA FOR CYLINDER FOLLOWING OGIVE IN SUB- 
SONIC AND TRANSONIC RANGE ARE INPUT TO PROGRAM or NO FIXED INPUT 
DATA FOR BO ATT AIL FOLLOWING LONG CYLINDER IN SUBSONIC AND TRANSONIC 
RANGE ARE INPUT TO PROGRAM is printed. 

For component 6 several estimated curves are used in the fixed input data. If these 
estimated curves play a major role in determining DA values for the vehicle, then the 
message THE D LAMBDA VALUES FOR THE FRUSTUM ARE ESTIMATED FOR THIS 
MACH NUMBER is printed. The various checks which the program performs on these 
estimated curves are described in detail in the section "Program Logic" found in the 
Programers' Section of this report. 

If a vehicle has a geometric parameter, or if a Mach number for which DA values 
are desired does not lie within the ranges given in table I, then errors may be introduced 
as a result of the extrapolation, and the results should be reviewed. 

If punched card output is obtained, the cards are punched exactly as the ju i coordi- 
nates and Da values are printed. The value of the ju coordinate will end in column 16, 
and the DA value will end in column 32. The format for the punched cards is 2E16.8. 

Two example configurations are shown in figure 5, and the computer printout for these 
examples is shown in table HI. The printout is shown in graphical form for each example 
in figure 6. 
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I 



f = 165 - 156-26? = 1-0 
e 8.735 


/^l\ = 8.735 

\ D o/l 11.513 


• 759 


(a) Example configuration 1. 



(b) Example configuration 2. 

Figure 5.- Two example configurations. All stations and diameters are in feet (meters). 
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TABLE III.- COMPUTER OUTPUT FOR EXAMPLE CONFIGURATIONS 1 AND 2 


EXAMPLE CONFIGURATION 1 

INPUT PARAMETERS 

NCMPS=5 NMACHN=2 NOXS=44 NDX=1 NC1=0 NC2=0 NC3=1 NFCYCM=1 NTHFRR=1 NPUNCH=0 
THNR=20 . 0 THBR=0. FAR=1. FNR=0. 

XR 4=0. XR5=0- XR7=25- DXh=0. DX5=0- DX7=8. XRlC=l6^. 


THFRR ( 1 ) 
12.500 


DDRR ( I ) 
.759 


COMPnO ( I ) 

DO (T) 

NXS(I 

3 

6.000 

3 

1 

8.735 

5 

5 

8,735 

8 

6 

8.735 

7 

7 

11.613 

21 


MACHs 1,000 


STATION 

1.75000000E+02 
1.7A0U0O00E+02 
1.73000000E+02 
1.71000000E+02 
1.69000000E+02 
1.67000000E+02 
1 .65000000E+02 
1 .65000000E+02 
1 .64000000E+02 
1 .63000000E+02 
1.62000000E+02 
1 .61000000E+02 
1 .60000000E+02 
1.58000000E+02 
1 .562b5000E+02 
1 . 5626500 OE + 02 
1 .55000000E+02 
1 .54000000E+02 
1 .53000000E+02 
1 .52000000E+02 
1.51000000E+02 
1.50000000E+02 
1 .50000000E+02 
1 . 48000000E + 02 
1.460000 00E+02 


dlamboa 

0 . 

2.20497778E+01 
1.110b35b6E+0l 
6.77071 547E + 00 
0.755031 1 3E + 0 0 
9.9 3 50719 6E +00 
8.21090000E+00 

1 .747O0OOOE+O0 

7.39400000E+00 
1 . 10916000E+01 
1.3850050OE+01 
1.56452500E+01 
1 .59067500E+01 
1.04436000E+01 
5. 15365000F. + 00 
2.44580000E+00 
4.521 16963E+00 
5.95017724E+00 
7 . 124434 1 2E + 0 0 
ft. 14288313E+00 
9.14230789E+00 
9.95874857E+00 
6.65323468E+00 
3.04494097E+00 
4.78672870E+00 
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TABLE III.- COMPUTER OUTPUT FOR EXAMPLE CONFIGURATIONS 1 AND 2 


1.44000000E+02 
1 .42000000E+02 
1 .40000000E + 02 
1.38NO0000F+02 
1*36000 OD 0 E + 0 2 
1.34000000E+02 
~T • TO 0 tro 0 0 0E+ 02 

' I.H5000000E+02 
1 .200D0000E + 02 
1.15000G00E+02 

" r.ioooooooE+o?' 

1 .OOOOOOOOE + 02 

7.50000000E+01 

5. 000 0000 0E + 0 1 “ 

2.50000000E+01 

2.50000000E+01 

1.500D000OE+01 

5.00000000E+00 

-5.00000000E+00 

MACH= 2.000 
STATION 

1 .75000000E+02 
1.74000000E+02 
1.73000000E+02 
1.71000000E+02 
1.69000000E+02 
1 .67000000E+02 
1.65000000E+02 
1.65000000E+02 
1 .64000000E+02 
1.63000000E+02 
1.62000000E+02 
1.61000000E+02 
1.60000000E+02 
1 .58000000E+02 
1 • 562650 00E+ 02 
1.56265000E+02 
1.55000000E+02 
1 .54000000E+02 
1.53000000E+02 
1.52000000E+02 
1 • 5 1 00 0 0 0 0 E + 0 2 
1 .50000000E+02 
1.50000000E+02 
1.48000000E+02 


4 • 51 005265E*00 
8.14148950E-01 
-1.64169983E+00 
-3.01453605E+00 
-3.51333 0 88E + 0 0 
-3.32138757£*00 
-1 •66516955E+00 
2. 294 77250 E- 02 
3.06695442E-01 
5.29463792E-01 
7 • 7 0 4 6 3 7 9 2 E - 0 1 
8.22357117E-01 
7.61584950E-01 
7.61584950E-01 
7.61584950E-01 
5.29200000E-01 
5 • 2 9 2 0 0 0 0 0 E - 0 1 
5.29200000E-01 
5.29200000E-01 


DLAM6DA 

0 . 

1 .44963333E+01 
8.61413048E+00 
8.38696873E+00 
1.11041449E+01 
1.38914817E+01 
1 • 61 59750 0E + 0 1 
6.98800000E+00 
6.98800000E+00 
6.9627O000E+00 
6.90065000E+00 
6.79945000E*00 
6.53620000E+00 
5.935 OOOOOE+OO 
5.24100000E+00 
8.38431544E+00 
9. 14661877E+00 
9.73594406E+00 
1 .03076740E+01 
1 • 090071 33E*0 1 
1.15083991E*0l 
1.21062634E+01 
2.38191 Q46E+00 
2.16431528E+00 


- Continued 
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TABLE III.- COMPUTER OUTPUT FOR EXAMPLE CONFIGURATIONS 1 AND 2 


1 .46000000E+02 
1.44000000E+Q2 
1.42000000E+02 
1.40000000E+02 
1.38000000E+02 
1.36000000E+02 
1.34000000E*02 
1 .300O0000E+02 
U25000000E + 02 
1 .20000000E+02 
1 ,15000000E+02 
1.10000000E*02 
1.00000000E+02 
7.50000000E+01 
5.00000000E+01 
2.50000000E+01 
2.50000000E+01 
1.50000000E+01 
5.00000000E+00 
-5.00000000E+00 


1.98250747E+00 
1.77213652E+00 
1 • 66351 143E + 00 
1.53182462E+00 
1.43029183E+00 
1.32372837E+00 
1.23023343E+00 
1.08457187E+00 
9.46492291E-01 
8.73043688E-01 
7.42800854E-01 
6.65644928E-01 
5.64452710E-01 
4.60654867E-01 
4.60654867E-01 
4.60654867E-01 
3.20093714E-01 
3.200937I4E-01 
3.20093714E-0I 
3.2n093714E-01 


- Continued 
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TABLE III.- COMPUTER OUTPUT FOR EXAMPLE CONFIGURATIONS 1 AND 2 - Continued 
EXAMPLE CONFIGURATION 2 

INPUT PARAMETERS 

NCMPS=4- NMACHN=2 NOXS=29 NDX=0 NC1=0 NC2=0 NC3=0 NFCYCM=1 NTHFRR=1 NPUNCH=0 

THNR=0. THBR=5* FAR=0. FNR=1.2 

XR^=0. XR5=0. xrt=o. DX^=0. DX5=0- dxt=o- XR1C=10. 

THTWrlT — — 

Q - - 


DORR (I) 

0 . 


COHPMOTI ) DO ( I ) >JXS ( I ) 

3 8.735 3 

2 10.000 6 

4 10.000 12 

8 10,000 8 


hlACHs 1.500 


STATION 

0 • 

1 .OOOOOOOOE+OO 
2. OOOOOOOOE+OO 
' 3 % OOOOOOOOE+O0 
4. OOOOOOOOE+OO 
6. OOOOOOOOE+OO 
8. OOOOOOOOE+OO 
1 .OOOOOOOOE+Ol 
1.20000000E+01 
1.40000000E+01 
1 .60000000E+01 

l.aoooooooE+oi 

2. OOOOOOOOE+Ol 
2.50000000E+01 
3.00000000E+01 
3 « 5 0 0 0 0 0 0 0 E ♦ 0 1 
4.00000000E+01 
5. OOOOOOOOE+Ol 
6. OOOOOOOOE+Ol 
6» OOOOOOOOE + O 1 
6,20000000E+01 
6.40000000E+01 
6.60000000E+01 
6.80000000E+01 
7» OOOOOOOOE+O 1 
7.50000000E+01 
8. OOOOOOOOE+Ol 


DLAM8DA 

0 . 

2,04935586E+01 
2 • 1 142 1 467E + 0 1 
1.90481776E+01 
1 .98268914E+01 
1.84077951E+01 
1 • 49709627E + 0 1 
8.14736272E+00 
5.51026129E+00 
4.83253871E+00 
4.20887419E+00 
3#6 392677 4E +00 
3.06966129E+00 
2.21375806E+00 
1.60112903E+00 
1.31708065E+00 
7.68080645E-01 
2.50500000E-01 
3.63709677E-02 
-7.50000000E-02 
-6.00719453E-01 
-1.39503351E+00 
-2,39868981E+00 
-2.95210336E+00 
**3» 1764650 2E +00 
-2.06523457E+00 
-1 #3489817 8E +00 
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TABLE III.- COMPUTER OUTPUT FOR EXAMPLE CONFIGURATIONS 1 AND 2 


MACHs 2.000 
STATION 

0 . 

1.00000000E+00 
2.00000000E+00 
3.00000000E+00 
4.00000000E+00 
6.00000000E+OQ 
8. OOOOOOOOE+OO 

1 . QOOOOOOOE+Ol 
1.20000000E+01 
1.40000000E+01 
1 .60000000E+01 
1.80000000E+01 

2. Q0000000E+01 
2.50000000E+01 
3.00000000E+01 
3.50000000E+01 
4. OOOOOOOOE+O 1 
5.00000000E+01 
6.00000000E+01 
6.00000000E+01 
6.20000000E+01 
6.40000000E+01 
6.60000000E+01 
6*8000000 0E +01 
7*000 OOOOOE+Ol 
7.50000000E+01 
8.00000000E+01 


DLAM8DA 

0 • 

1 • 781 56951E+0 1 
1 .88320817E+01 
1 .74064891E+01 
1.76526230E+01 
1.48578013E+01 
1.14253912E+01 
6.55813882E+00 
8.27437419E+00 
4.66842S81E+00 
4*1091 1613E+00 
3.59644516E+00 
3.08377419E+00 
2.27383871E+00 
1 *68741 935E+00 
1.38038710E+00 
8.81387097E-01 
3.88000000E-01 
1.62580645E-01 
-3.00000000E-01 
-6.92395514E-01 
-1.06952569E+00 
-1.38506436E+00 
-1 .67063679E+00 
-1.92444277E+00 
-1 .43828836E+00 
-7.96386833E-01 


- Concluded 
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PROGRAMERS' SECTION 


PROGRAM LOGIC 


1 dCjj 

In referring to reference 1, note that the fixed input data are values of S — - — —(= X). 

D dx 

Hence, after the interpolations on X are performed, the program multiplies the results 

dCN 

by D to obtain the final form S S.(= DX). This process is applied to every case 

dx 

input x station for each component comprising the vehicle. 


For best results the pro- 
gram should be used for vehi- 
cles which have geometric 
parameters within the ranges 
given in table I, and the case 
input Mach numbers should, if 
possible, also be within the 
ranges given in table I. If 
certain parameters or Mach 
numbers are outside these 
ranges, then the program may 
still be used to obtain DX 
distributions with the under- 
standing that sizable errors 
may be introduced as a result 
of the extrapolation. In order 
to decrease the size of these 
errors, X values for limiting 
values of certain parameters have been incorporated as a part of the program. These X 
values in general do not represent experimental results, and anytime these values are 
used in an interpolation, the final results should be reviewed as to their reliability. The 
geometric significance of the limiting parameter values are shown in figure 7. The 
description of the X values which correspond to the limiting parameter values are given 
in table IV. Some of the limiting parameters and associated X values are inputs to the 
program and some are calculated in the program. 

In obtaining the X values for an ogive with f n = 0.5 one can use with some modi- 
fication the equation for the X values for a sphere given by equation (1). The value of 
A can be obtained for the desired Mach number by interpolating the proper A values 




Figure 7.- Geometric significance of parameter limitations for 
components 6 and 7. 
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TABLE IV.- PARAMETER LIMITATIONS AND X VALUES USED FOR THESE LIMITATIONS 


Component 

Parameter 

limitation 

X values used for the limitations 

1 

e n = °° 

X is taken to be a small number, usually 0.01 or 0.001 for each x/L 
station and Mach number * 

2 

f n = 0.5 

When f n = 0.5, the ogive becomes a sphere; hence the X values were 
obtained from sphere values * 

3 

None required 


4 

f n = 0.5 

No fixed input data in program for cylinder following sphere 

5 

D 

11 

O 

o 

X = 0.01 for all x/L stations and Mach numbers except Mach 15 in 
which case X = 0.116 for all x/L stations* 

6 

6f = 0° 

X values are obtained from the X values for a cylinder following a 
cone beginning at station XI as shown in figure 7t 


f a = o 

X values are obtained from X values for a cone* 

7 

o 

? 

(1 

o 

X values are obtained from X values for a cylinder following a conet 


D lPo = 1 or 

X values are obtained from a cylinder following a cone beginning at 


6f = 0° 

station X2 as shown in figure 7^ 

8 

e b = o° 

X = 0.0 for all x/D 0 stations and Mach numbers* 


*These X values are input to the program, 
t These X values are calculated in the program. 


given in the fixed data listing for sphere data. The modified equation is 


where x/L assumes the values given in the fixed data listing for sharp tangent ogives. 

The A values for hypersonic Mach numbers were calculated by means of the modi- 
fied Newtonian theory as presented in reference 3. These values are input to the program 
as Mach 15 data and serve as an upper boundary on Mach number. For any case input 
Mach number greater than or equal to 15, the program will use these values and no inter- 
polation on Mach number will be performed. 

The program can accommodate five variations from the list of eight basic components 
shown in table I. These variations and associated parameters are shown in figure 4. In 
obtaining the DA distributions for variations I, n, and ID; the program considers each 
variation as being a single cylinder, calculates the A values for this cylinder — it should 
be recalled that the A values are independent of body diameter - and multiples these 
values by the appropriate diameter. The DA results obtained for variation IV are valid 
only when f a = 4.0. Again the A values associated with these variations do not repre- 
sent experimental data, and the DA distributions obtained for any of these variations 
should be reviewed. 
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A simplified flow chart for the program is shown in figure 8. Detailed flow charts 
for the individual components are shown in figures 9 to 11. 

The process which the program uses to obtain the DX distributions for compo- 
nents 1, 6, and 7 is shown in figures 9 to 11, respectively. In general the process consists 
of linearly interpolating the proper fixed input data values. If no interpolation is required 
at some of the case input parameters, then the program will omit some of the steps shown 
in figures 9 to 11. If insufficient fixed input data is present to permit an interpolation, 
then the program will extrapolate. 

The X values for components 2, 4, 5 and 8 are functions of three parameters - 
Mach number, x/D 0 (x/L for component 2), and the appropriate other parameter 
(f n for components 2 and 4, e n for component 5, and 0^ for component 8). Conse- 
quently, the steps used to obtain the DX distributions for each of these components are 
similar to those used for component 1. 

The DX values for component 3 are obtained by use of the following closed-form 
equation: 

DX = 8*D 0 A(0.5 - £)(l - £)(£) W 

which is derived from equation (A16) of reference 1. For a given case input Mach num- 
ber, the program will interpolate on the fixed input A values to find a value to use in 
equation (1), will convert the case input x stations into x/D 0 values, and will substi- 
tute these values into equation (1). 
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Figure 9.- Flow chart for component 1. 
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COMPUTER PROGRAM 


The computer program consists of a main program, three main subroutines, and the 
fixed input data. The name and purpose of each subroutine is described as follows: 


Subroutine 


Purpose 


MACHXD Converts case input ji coordinates into x/D 0 

or x/L stations, determines if interpolation 
at case input Mach number is required and if 
so finds two fixed input Mach numbers to use 
in interpolation 

INTP Determines DX values for components 1, 2, 4, 

5, and 8 

DISCOT Performs all interpolations for all components; 

sometimes this subroutine is used to perform 
simultaneous double interpolations at two inde- 
pendent variables, and sometimes it is used to 
interpolate at only one independent variable 


The DISCOT subroutine includes three subroutines, UNS, DISSER, and LAGRAN, 
which are never called for in the main program but are ah inherent part of the DISCOT 
subroutine. 


The purpose of the main program is to read in the fixed input data, read in the case 
input data, convert double or triple subscripted variables into single subscripted variables, 
and to obtain DX values for components 3, 6, and 7. 

The programing is FORTRAN IV. 


METHOD OF INPUTTING FIXED INPUT DATA 


A complete listing of all the fixed input data which are input to the program is 
enclosed. These data include the following: (1) code numbers - these numbers, with one 
exception, tell the program "how many" parameter and x values are being input. The 
exception is the set of code numbers which tell the program whether or not the frustum 
X values are experimental or estimated; (2) the values of the parameters; and (3) the 
X values. 

The order in which these data appear in the data listing is the order in which they are 
input to the program, and this order must always be preserved. The cards which give 
the FORTRAN names of the variables are listed solely for that purpose and are not a part 
of the data. The definitions of these FORTRAN names can be found in table V. 

The method of inputting the data is best described by considering each component 
separately. 
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TABLE V.- FORTRAN NAMES AND DEFINITIONS FOR FIXED INPUT DATA 


FORTRAN name 

Definition 

FORTRAN names for fixed input data for component 1 

NMACH(l) 

Number of Mach numbers 

NXL1 

Number of x/L values 

NTHNl(I) 

Number of (9 n values for Ith Mach number; I = 1 to NMACH(l) 

NLAMl(I) 

Number of A values for Ith Mach number; I = 1 to NMACH(l) 

MACH(1,I) 

Values of the Mach numbers; I = 1 to NMACH(l) 

XL1(J) 

x/L values; J = 1 to NXL1 

THN1(I,K) 

0 n values at the Ith Mach number; I = 1 to NMACH(l), K = 1 to NTHNl(I) 

LAM1(I,L) 

A values at the Ith Mach number; I = 1 to NMACH(l), L = 1 to NLAMl(I) 


FORTRAN names for fixed input data for component 2 

NMACH(2) 

Number of Mach numbers 

NXL2 

Number of x/L values 

NFN2(I) 

Number of f n values for Ith Mach number; I = 1 to NMACH(2) 

NLAM2(I) 

Number of A values for Ith Mach number; I = 1 to NMACH(2) 

MACH(2,I) 

Values of the Mach numbers; I = 1 to NMACH(2) 

XL2(J) 

x/L values; J = 1 to NXL2 

FN2(I,K) 

f n values at the Ith Mach number; I = 1 to NMACH(2), K = 1 to NFN2(I) 

LAM2(I,L) 

A values at the Ith Mach number; I = 1 to NMACH(2), L = 1 to NLAM2(I) 

FORTRAN names for fixed input data for component 3 

NMACH(3) 

Number of Mach numbers 

MACH(3,I) 

Values of the Mach numbers; I = 1 to NMACH(3) 

A(I) 

A values; I = 1 to NMACH(3) 

FORTRAN names for fixed input data for component 4 

NMACH(4) 

Number of Mach numbers 

NXD4 

Number of x/D 0 values 

NFN4(I) 

Number of f n values for Ith Mach number; I = 1 to NMACH(4) 

NLAM4(I) 

Number of A values for Ith Mach number; I = 1 to NMACH(4) 

MACH(4,I) 

Values of the Mach numbers; I = 1 to NMACH(4) 

XD4(J) 

x/D 0 values; J = 1 to NXD4 

FN4(I,K) 

f n values at Ith Mach number; I = 1 to NMACH(4), K = 1 to NFN4(I) 

LAM4(I,L) 

A values at Ith Mach number; I = 1 to NMACH(4), L = 1 to NLAM4(I) 
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TABLE V.- FORTRAN NAMES AND DEFINITIONS FOR FIXED INPUT DATA - Continued 


FORTRAN name 

Definition 


FORTRAN names for fixed input data for component 5 

NMACH(5) 

Number of Mach numbers 

NXD5 

Number of x/D 0 values 

NTHN5(I) 

Number of 9 n values for Ith Mach number; I = 1 to NMACH(5) 

NLAM5(I) 

Number of X values for Ith Mach number; I = 1 to NMACH(5) 

MACH(5,I) 

Values of the Mach numbers; I = 1 to NMACH(5) 

XD5(J) 

x/D 0 values; J = 1 to NXD5 

THN5(I,K) 

d n values at Ith Mach number; I = 1 to NMACH(5), K = 1 to NTHN5(I) 

LAM5(I,L) 

X values at Ith Mach number; I = 1 to NMACH(5), L = 1 to NLAM5(I) 


FORTRAN names for fixed input data for component 6 

NMACH(6) 

Number of Mach numbers 

NXD6 

Number of x/D 0 values 

NTHN6(I) 

Number of Q n values for Ith Mach number; I = 1 to NMACH(6) 

NTHF6(I,J) 

Number of 9f values at Ith Mach number and Jth 0 n ; 1 = 1 to NMACH(6), 

J = 1 to NTHN6(I) 

NFA6(I,J,K) 

Number of fa values at Ith Mach number, Jth 0 n and Kth 0f ; 
I = 1 to NMACH(6), J = 1 to NTHN6(I), K = 1 to NTHF6(I,J) 

NLAM6(I,J,K) 

Number of X values at Ith Mach number, Jth 0 n . and Kth Of, 
I = 1 to NMACH(6), J = 1 to NTHN6(I), K = 1 to NTHF6(I,J) 

NTOTL6(I,J) 

Total number of X values at Ith Mach number and Jth 0 n ; 
I = 1 to NMACH(6), J = 1 to NTHN6(I) 

NTOTFA(I,J) 

Total number of f a values at Ith Mach number and Jth 0 n ; 
I = 1 to NMACH(6), J = 1 to NTHN6(I). 

MACH(6,I) 

Values of the Mach numbers; I = 1 to NMACH(6) 

XD6(L) 

x/D 0 values; L = 1 to NXD6 

THN6(I,J) 

e n values at Ith Mach number; I = 1 to NMACH(6), J = 1 to NTHN6(I) 

THF6(I,J,K) 

6* values at Ith Mach number and Jth 0 n ; 1 = 1 to NMACH(6), 

J = 1 to NTHN6(I), K = 1 to NTHF6(I,J) 

FA6(I,J,M) 

f a values at Ith Mach number and Jth 0 n ; 1 = 1 to NMACH(6), 

J = 1 to NTHN6(I) , M = 1 to NTOTFA(I,J) 

K1(I,J,M) 

Code numbers denoting experimental or estimated data at Ith Mach 
number and Jth Q n \ I = 1 to NMACH(6), J = 1 to NTHN6(I), 

M = 1 to NTOTFA(I,J) 

LAM6(I,J,N) 

X values at Ith Mach number and Jth 0 n ; 1 = 1 to NMACH(6), 

J = 1 to NTHN6(I), N = 1 to NTOTL6(I,J) 
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TABLE V.- FORTRAN NAMES AND DEFINITIONS FOR FIXED INPUT DATA - Concluded 


FORTRAN name 

Definition 


FORTRAN names for fixed input data for component 7 

NMACH(7) 

Number of Mach numbers 

NXD7 

Number of x/D 0 values 

NDDR7(I) 

Number of Bi/D q values for Ith Mach number; I = 1 to NMACH(7) 

NTHF7(I,J) 

Number of Of values for Ith Mach number and Jth Dj Id 0 ; 
I = 1 to NMACH(7), J = 1 to NDDR7(I) 

NLAM7(I,J) 

Number of X values at Ith Mach and Jth Di/D 0 ; I = 1 to NMACH(7), 
J = 1 to NDDR7(I) 1 

MACH(7,I) 

Values of the Mach numbers; I = 1 to NMACH(7) 

XD7(K) 

x/D 0 values; K = 1 to NXD7 

DDR7(I,J) 

Df/D o values at Ith Mach number; I = 1 to NMACH(7), 
J = 1 to NDDR7(I) 

THF7(I,J,L) 

Of values at Ith Mach number and Jth Di/D n ; I = 1 to NMACH(7), 
J = 1 to NDDR7(I), L = 1 to NTHF7(I,Jr 

LAM7(I,J,M) 

X values at Ith Mach number and Jth Di /D 0 ; I = 1 to NMACH(7), 
J = 1 to NDDR7(I), M = 1 to NLAM7(I,J) 

FORTRAN names for fixed input data for component 8 | 

NMACH(8) 

Number of Mach numbers 

NXD8 

Number of x/D 0 values 

NTHB8(I) 

Number of 0 b values at Ith Mach number; I = 1 to NMACH(8) 

NLAM8(I) 

Number of X values at Ith Mach number; I = 1 to NMACH(8) 

MACH(8,I) 

Values of the Mach numbers; I = 1 to NMACH(8) 

XD8(J) 

x/D 0 values; J = 1 to NXD8 

THB8(I,K) 

0 b values at Ith Mach number; I = 1 to NMACH(8), K = 1 to NTHB8(I) 

LAM8(I,L) 

X values at Ith Mach number; I = 1 to NMACH(8), L = 1 to NLAM8(I) 


Component 1 

Corresponding to each Mach number, there is associated a X value array which is 
a function of the 0 n values for which data are available for this Mach number and also 
a function of a fixed set of x/L values. These X values are input columnwise as 
denoted by the arrows in table VI. The complete set of X values for all Mach numbers 
are input as a two-dimensional array where the first subscript refers to Mach number, 
and the second refers to the position in the associated X value array. 

The fixed set of x/L values, the Mach numbers, and the corresponding sets of 
9 n values are all input in ascending order. 

An example of the X value array form for a given Mach number is included in 
table VI. (Table VI is found at the end of the Programers' Section.) 
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Components 2, 4, 5, and 8 


The method of inputting the storage data for each of these components is similar to 
the method used for component 1. An example of the X value array form for a given 
Mach number for each of these components is included in table VI. 

Component 3 

The DA values for this component are determined by the closed-form equation 
(eq. (1)); hence, no X values are input. The only parameters input are Mach number 
and corresponding coefficients A(I). The Mach numbers are input in ascending order. 

Component 6 

Corresponding to each Mach number, 0 n , and Of combination, there is associated 
a X value array which is a function of the f a values for which data are available for 
the given Mach number, 0 n , and Of combination and also a function of a fixed set of 
x/D 0 values. Hence, for a given Mach number and Q n combination, there exists as 
many X value arrays as there are Of values which correspond to the Mach number and 
0 n value. These individual arrays are combined into a single array, which now - for a 
given Mach number and 0 n value - is a function of f a , Of, and x/D 0 . These X 
values are input columnwise as denoted by the arrows in table VI. The complete set of 
X values for all Mach numbers and 0 n values are input as three-dimensional variables 
where the first subscript refers to Mach number, the second refers to the 0 n value, 
and the third refers to the position in the associated X value array. 

The Mach numbers and corresponding sets of 0 n values are input in ascending 
order. The sets of Of values which correspond to each Mach number and 0 n combina- 
tion are also input in ascending order. The sets of f a values which correspond to each 
Mach number, 0 n , and Of combination are input in ascending order. The fixed set of 
x/D 0 values are input in ascending order. 

If for a given Mach number, 0 n , Of, and f a combination; the associated X value 
curve is estimated, then the value of the fixed input variable K1 is 0. If the associated 
X value curve is experimental, then the value of K1 is 1. 

An example of the array form for a given Mach number and 0 n value is included 
in table VI. 


Component 7 

Corresponding to each Mach number and Di j D 0 combination, there exists an 
array of X values which is a function of the Of values for which data are available for 
this Mach number and /D q combination and also a function of a fixed set of x/D 0 
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values. These X values are input columnwise as denoted by the arrows in table VI. The 
complete set of X values for all Mach numbers and Di/D 0 combinations are input as 
three-dimensional variables where the first subscript refers to Mach number, the second 
refers to the Dj/D 0 value, and the third refers to the position in the associated X 
value array. 

The fixed set of x/D 0 values, the Mach numbers, and the corresponding sets of 
Dj/Dq values are all input in ascending order. The df values which correspond to each 
Mach number and Dj/Dq combination are also input in ascending order. 

An example of the array form for a given Mach number and Dj/D 0 combination 
is included in table VI. 
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TABLE VI.- EXAMPLE OF FORM OF X-VALUE ARRAY FOR INPUTTING FIXED 


Component 


1 


2 


4 


INPUT DATA FOR EACH COMPONENT 


XL1(J) 
THN(1,K) 

0 

10 

15 

25 

33 

40 

XL2(J) 


Example of form of X-value array 
Mach(l,l) = .7 


0.0 

LAM1(1,1) = .00 
LAM1(1,2) = .00 


LAM1(1,6) = .00 


.2 

LAM1(1,7) = .01 


.4 .6 .8 


1.0 


LAM1(1,31) = .01 


LAM1(1,36) = .35 


Mach(2,l) = .8 


FN2(l,ig''\^ 

0.0 


.05 .10 - 

1.00 

.50 

LAM2(1,1) = 0.00 

LAM2(1,5) = 2.58 

LAM2(1,8 

1) = 0.00 

1.39 

LAM2(1,2) = 4.15 





3.00 

LAM2(1,3) = 1.95 



\ 

/ 

6.00 

LAM2(1,4) = 1.25 

\ 

t 

LAM2(1,84) = 0.01 

The complete set of XL2(J) for each array for this component is . 

,00, .05, .10, 

.15, .20, 

.25, .30, .35, 

.40, .45, .50, .55, .60 ; 

, .65, .70, 

.75, .80, .85, .90, .95 

, 1.00. 




Mach (4,1) = 3.0 



^^\XD4(J) 






FN4(1,K)^^ 

0.0 


.5 1.0 - 

6.5 

3 

LAM4(1,1) = .405 

LAM4(1,4) = .320 

LAM4(1,40) = .025 

5 

LAM4(1,2) = .270 



4 


7 

LAM4(1,3) = .178 

\ 


LAM4(1,42) = .013 


The complete set of XD4(J) for each array for this component is 0.0, .5, 1.0, 1.5, 2.0, 
2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 5.5, 6.0, 6.5. 

Mach(5,2) = .8 



0.0 

LAM5(2,1) = .01 
LAM5(2,2) = .63 


6.5 


LAM 5(2, 8) = .01 


LAM5(2,7) = -2.47 


LAM5(2,162) = .01 


LAM5(2,168) = .03 


The complete set of XD5(J) for each array for this component is .0, .1, .2, .3, .4, .5, .6, 
.7, .8, .9, 1.0, 1.25, 1.5, 1.75, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 5.5, 6.0, 6.5. 
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TABLE VI.- EXAMPLE OF FORM OF X- VALUE ARRAY FOR INPUTTING FIXED 
INPUT DATA FOR EACH COMPONENT - Concluded 


Component 


Example of form of X- value array 


XD6(L) 
FA6(7,1,mT 


Mach(6,7) = 4.04, THN6(7,1) = 15°, 
THF6(7,1,1) = 5° 


0.0 


0.2 


0.4 


6.5 


LAM6(7,1,1) = .52 
LAM6(7,1,2) = .74 
LAM6(7,1,3) = .30 


LAM6(7,1,4) = .52 


THF6(7,1,2) = 20° 


LAM6(7,1,52) 

LAM6(7,1,54) 


1 [ LAM6(7, 1,55) = 2.05 LAM6(7,1,57) = 1.90 LAM6(7,1,89) 

4 | LAM6(7,1,56) = 1.80 4- LAM6(7,1,90) 

The complete set of XD6(L) for each array for this component is 0.0, 0.2, 0.4, 0.6, 0.8, 
1.25, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 5.5, 6.0, 6.5. 


Mach(7,l) = .80, DDR7(1,1) = .70 


XD7(K) 

thf7(i,i4lT 


0.0 


0.1 


0.2 


6.5 

LAM7(1,1,93) 


5 LAM7(1,1,1) = -.370 LAM7(1,1,5) = .070 

10 LAM7(1,1,2) = -1.140 

15 i , 

30 LAM7(1,1,4) = -1.500 i LAM7(1,1,96) 

The complete set of XD7(K) values for each array for this component is 0.0, .1, .2, .3, A 
.6, .7, .8, .9, 1.0, 1.25, 1.5, 1.75, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 5.5, 6.0, 6.5. 


XD8(J) 


THB8(1,K) 


0 

4 

8 

12 

16 


Mach(8,l) = 1.50 

0.0 0.1 


.2 - - 


2.5 


LAM8(1,1) = .000 LAM8(1,6) = .000 

LAM8(1,2) = .000 


LAM8(1,56) = 


LAM8(1,60) = 


LAM8(1,5) = - .060 

The complete set of XD8(J) values for each array for this component is 0.0, .1, .2, .4, .6 
1.0, 1.25, 1.5, 1.75, 2.0, 2.5. 


...i 

= .52 
= .45 

= 2.30 
= 1.90 

1 . 0 , 

= .06 

= .095 
1, .5, 

.000 

-.060 

, . 8 , 
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CONCLUDING REMARKS 


The results obtained from this program yield adequate characteristics for most 
launch- vehicle configurations at angles of attack less than 5°. Because of the lack of 
sufficient experimental data which are input to the program, certain discrepancies exist 
between experimental and empirical results. For components which have geometric 
parameters within the ranges of the fixed input data, these discrepancies are small. For 
components which have geometric parameters outside these ranges, the discrepancies 
may be larger, and the empirical results should be interpreted only as representing the 
trend of the characteristics curves. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., August 21, 1967, 

124-08-05-25-23. 
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APPENDIX A 


PROGRAM LISTING 



£> REAL MACH , MR ,MR6, MR7, Ml ,LCN , MACHN, LAM 1 , LAM2 , LAM4 , LAM5 »LAM6 , 

05 XLAM7 * LAN'Si LAM 

i‘ntege^”c6mpno“ 

EXTERNAL DISCOT 

COMMON NMACH (8)_,NTHN1 ( 1 4 ) , NLAM l_( 1 4J , MACH ( 8 , 1_6 )_, THN1 ( 1 4 , 1 7]_,_XL_1 ( 6 L’ 
1 L AM 1 (14,102) ,NL AV2 ( 7 )' ,'NFN2~( 7*7’, XL2< Z \ ) , FN2 (7 , 5 ) ,L A M2 T'i , \ 05 ) , a"< 127, " 
2NLAM4 ( ^ ) ,njFN4 ( « ) , XD4 ( 14 ) ,FN4 (5,3) ,LAM4 (“5,42 ) , NLAM5 (16) , N THN5 (16), 
3XD5 (24 ) * THN5 (16*9) ,LAM 5(16,216) , NT HN6 ( 3 ) , NTHF6 (8,3), 

4NF A6( 8 , 3,4 ) ,'NL AM6 "( 8 ,' 3",'4) , NT0TL6( tT,''3 ) , N 1 0 T Fm”( ^3V,ThKT 6( bTSTI 
5XD6 (18) , THF 6 ( 8 , 3 , 4 ) , L A M 6 (6,3,288) ,FA6 <8 ,3 , 1 6 ) ,.<1 (8,3, 16) , 

6NDDR7 ( 8 ) ,NTHF7 (8,3) ,NLAM7(8 , 3 ) , DDR7 (8,3) , TMF7 (8, 3, 4 ) ,LAM7( 8 ,3,96 ) , 

7X D 7' ( 24 )' , N T HB6 (' S ) V'THc 8( S ib ) 7X08 ( TZ T , L AM 8 < 8%' 60 1 7iNL AMST 8 ) 

COMMON HEAD ( 1 1 ) ,XX( 100) , X ( 25 ) , LAM ( 228 ) , NXS ( 14 ) »X3(350 ) , DO ( 14) , 

1 A 1 (2 , 1 00 ) , VA ( 1 7 ) , TN (2 ) , MACHN (15) 

COMMON A2 ( 2 , 100 ) , A 3 (2,1 00), M 1 " (*2 ) ~, A I l" (27TA4“<T00 ) » AA1 '< 1 CO )'♦ 

1 AN3 (5 ) ,XX6 ( 1 00 ) ,FA (4 ) ,i< 1 i (4 ) , THF (4 ) , A22 (2 ) ,A33 (2 ) ,XX7 ( 1 00 ) , AA2 (2 ) , 
2C2 (2 ) , A 1 2 (2 ) , CQMPNO ( l a ) , A44 (2,100), THFRR ( 5 ) , DDRR (5 ) 

C READ STATEmENTS’ eor' E l XED Fn PUT 5 A T A~ 'f 0 R~ SHARP "CONES 

READ ( 5 , 9502 ) NMACH ( 1 ) 

READ ( 5 ♦ 9502 ) NXL 1 

NM 1 =NM A CH ( 1 ) ' " ' ' " ' ‘ 

READ (5, 9502 ) (NTHN1 (I ) , I =1 , N M 1 ) 

READ ( 5 , 9502 ) (NLAM 1 ( I ) , I = 1 , NM 1 ) 

READ ( 5,1 60 J ) (MACHTlTl ) ,T=T,KFTl ) 

READ (5, 1600) ( XL 1 ( J) , J = 1 ,NXL1 ) 

DO 1 501 1=1, NM i 

NTN'ZnTHNM II ’ 

1 5 r 1 READ (5, 1 5 n 0 ) ( THN 1 (I , K ) , <= 1 ,NTN ) 

001 50?I = 1 , NM l 

NL = N|_ A v 1 ( I j " " ' 

1502 READ (5 , 1 6^0 ) (LAV 1 ( t , l ) ,L=1 , NL ) 

C REAP STATEMENT* ^OR FIXED INPUT DATA FOR SHARP TANGENT OGIVES 
READ (5', 9 50 2") nMT^hT 2") 

READ ( 5 , 9502 ) NXL 2 

NM2 = NVACH (2 ) 

READ ( 5 , 95 bZ) (N e ' N2~( T)7i=T,N7 2 ') 

READ(5,9502 ) (NLA - v 2 ( I ) , I =1 ,NV? ) 

READ(5, 1600) (MACH(2, i ) , I = r,NV2) 

READ (5, 1600 ) (XL2 ( Jl7 J = r<T\XL2 ') " " 
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DO 1 503 I = 1 ,NV2 
NFN=NFN? ( I ) 

15~2 READ (5 « 1600 ) (-M2 ff ,<)»< = 1 * \-N ) '* ' 

DC1504I =1 ,NV!2 
ML =N'L A v2 ( I ) 

15 r 4 READ ( 5 * 1500 ML'Avp (1 V|_ ) ♦ L=1'»NL') ~~ ““““ 

C READ STATEMENTS FOR FIXED i-\PJT DATA FOR SPHERES' 

READ (E , 95**2 INiVACH (3 ) 

\i f /3 = \iv ACM ( 3 ) 

READ(5» 1600) (MACH (3« I) » ! = 1 iW3 ) 

REAC (5, 1 600 ) ( A ( I ) » ! = 1 , NfO ) 

C ’READ SYiTf *VtN Tr=OP"r IXcD''"' INPUT "DAT A T3R~ CYLINDERS FOLLOWING TANGENT 
C OOIVE c 

READ (5, 9 C 02 )NMACH (4 ) 

R - A 0 ( 5 Y 9 E'hF j'\ x O 4 
\! v 4 = M v A C H ( 4 ) 

READ'(h,'95h2 ) (N-N4 ( I ) ♦ I = 1 ,N:VA ; 

r’eadT =795*^*) Tnla t: a"(T r i" i”=TTn xaT 

READ ( 5 * 1 60C ) (MACn ( 4 , I ) , I = 1 i N,M4 ) 

READ (5il 60U ) C XD4 ( J ) « J= 1 « NXD4 ) 

CO 1 507T'=’iTnm4 * 

NFN = NFN4 ( I ) 

1 507’ READ (5, 1600) (FN4 ( I « k ) , <=1 iNFN ) 

DO 1 5 'M3 I = I , JM4“ ‘ 

NL=N|_AM4 ( I ) 

1508 READ ( 5 i 1600) (LAIX4 ( I «L) iL=l iNL ) 

C READ S TA T E i^ENT S“F 0 R " ' FT X E” I N P ■ JT D A TA" — FO R ' CYLINDERS'" FOLLOW I NG CONES 
READ ( 5 i 95H2 ) MmACH ( 5 ) 

READ (5» 95^2 )NXD5 ' ' 

VM 5 = NI m A CH( 5 ') 

READ (5, 9502 ) (NTHN5 ( I ) « I = 1 »NM5 ) 

READ(5<9502) (NLAM5 ( I ) *1=1 » NM5) 

READY 5 , 1 600 ) ( MACH 7*5 • I f» 1 = 1 » NM5 ) ~ 

READ (5, 1 600 ) ( XD5 ( J ) * J = 1 ,NXD5 )_ 

DO 15051=1 » NM5 ----- 

NT N = N THNSTi') 

1505 READ ( 5 » 1600) (THN5( I t<) » < = 1 « NTN ) 

DO l 506 I = 1 «NM5 ” ’ 

*>• nL=NLAV5(T’)’ " " ” ~ ’ ” 
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00 


1 506 READ (=5, 160 0 ) (LAME ( [ ,L ) » L= 1 t NL ) 

C READ STATEMENTS FOR FIXED INPUT DATA FOR FRUSTUMS 'following CONE" 
C 'CvLTnOERS ’ 

READ ( 5 , 9502 ) NVACH (6 ) 

READ ( 5 * 9502 ) NXQ6 

NM6=NMACH ( 6 ) 

READ ( 5» 9502 ) (NTHN6 ( I ) » I = 1 »NM6 ) 

DO 1 509 I = 1 ,NM6 
N t N = NT HN 6 ( I ') 

15~9 READ (5, 95^2 ) ( N T H F 6 ( I ♦ J ) * J= 1 * NTN ) 

D01 ^1 ° T = 1 ,NM6 ' 

NT N ="nT H N6TI ’ ) 

DO 1 C 1 ^J = l ,MTM 
MT^ = \iTHF6 ( I , J ) 

1 5 1 r ' 'Re a DTo T^dST: VY \"F A 6TTTJ ;<y VxT= 1 «NTF ) 

DO 1 5 1 1 1 = 1 * NMS 
N T N = Ki’T HN6 ( I ) 

D 5T5TT J=T TnTn 

NTF = NTHF6 ( I 1 J ) 

1511 READ ( 5 » 9502 ) (NLAiV6( I , J » < ) U<= 1 « NTF ) 

D01 51 21 = 1 ,NM6 
NTN=NTHN6 < I ) 

1512, REA D ( pj 950 2 ) ( N I T 0 T i_ 6_( 1 _» J ) i J= 1 « NT N) 

' DO l"513T= r."jM6 - - - - - - - 

NTN=NTHN6 ( I ) 

1513 READ ( 5 , 9502 ) (NTOTFA ( I , J ) , J= 1 ,NTN ) 

READ" ( 5 1 1 60": r n MACH ( 6 i f) , I = i ♦ NXI6~ ) 

READ ( 5 i 1 600 ) ( X D 6 ('L ) *L = 1 <NX06 ) 

D01 5 1 4 I = 1 ,NM6 
N'T N=NTHN6 ( 1 ) 

1514 READ (5*1 600 ) ( THN6 ( I i J ) * J = 1 i NTN ) 

DO 15151=1 ,NM6 

NTN=NTHN6(i) - — 

D0151 5 J = 1 i NTN 
NTF = \'THF6 ( I i J ) 

1515 READ ( 5 « 1600) (ThF 6 ( I » J i < ) i < = 1 « NT F ) 

DO 1 526 I = 1 «NM6 

NTN=NTHN6 ( I ) 

DO 1 526 J= 1 i NTN 
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NT OTF = NT CT~ A C I J ) 

1526 R E A D ( 5 i 1600) < r A6< i « J » \A ) »M=1 * N T O T F ) 

COl 5 1 6 I = 1 ,NW6 ' ’ 

NTN=NTHN6 ( I ) 

001 51 6J= 1 ,\TN 
NT0TF=NT0TFA { I , j ) 

15 16 READ ( 5 ♦ 9502 ) << 1 < I « J i v. ) , 1 « NT OTF ) 

DO 1 5? 7 ! = 1 , NMf 

NTN = mT u \ | 6 ( I ) 

DO 1 52 7 J = 1 , NT ’ i 
NT0TL = NT0TL6 ( I « J ) 

1527 READ (5, 1600) (LAV6( i iJiN) »N=1 ,\|TCTL‘> ” 

C READ STATEMENTS FOR FIXED INPUT DATA FOR CYLINDERS FOLLOWING FRUSTUMS 

READ (5, 9502 ) NMACH (7 ) 

READ ( 5 « 9502 ) NXD7 " ' 

N M 7 = M mi A C H ( 7 ) 

READ(5»95°2 ) ( NDDR 7 ( I ) i I =1 »NM7) 

C0152ST=1 .NM - ’ ~ 

ND = NDDR7 ( I ) 

1528 READ (5, 9502 ) (NTHF7( I , J ) *J=1 «ND) 

DC 1517 1 = 1 « NM?" " 

N0=ND0R7(I) 

1517 R-AD(5»9FD2 ) (NLAM7( I « J ) » J=1 «ND) 

READ (5, 1 6 0 J ) (MACH (7,1) , 1=1 t N M 7 ) ■ - - ■ 

READ ( 5 « 1600) (XD7(< ) *K = 1 .NXD7) 

CO 1 5 1 8 I = 1 ,NM7 

ND = NDDR7(I) ' * 

1518 READ (5 , 1 600 ) (DDR7 ( I « J ) « J= 1 « ND ) 

CO 1 5.29 1=1 « NM 7 

ND = NDDR7(I) ‘ '' 

CO 1 529 J= 1 < ND 
NTF = NTHF7 ( I « J ) 

15 29 R c A0 (5 , 160 0 ) (THF7 ( I , J , L ) < L = 1 , NTF ) 

0015191 = 1 * N V 7 — ' 

ND=NDDR7 ( I ) 

C 015 1 9 J =T * NO 
NL=NLAM7 ( I » J ) 

1519 READ (5, 1600 ) (LAM7 ( I « J *M ) «Vi=l , NL ) 
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U1 

o 


"C READ STXf £m?NtT“F 0 R ~T7^1jT ~ DATA' ~F0 R ‘ ~3Q A TT^TlITf OLi_ OWING" UO N(T 

C CYLINDERS 

READ t 5 » 9502 ) NMACH (8 ) 

"Read 1579’5o'2Tnx'd'8 

MM8=M^ACH(8) 

Read (5, 9502) (NTHee ( i ) » i = i » nmq ) 

Re ad ( 5%" 9eo? )'Tnl a ms" '( i ) 7 r=T7NM8l 

READ ( 5 » 1600) ( VACH ( 8. I ) , 1 = 1 »NM8 ) 

RE AD ( 5, 1 6^0 ) (XD8 ( J ) * J= 1 .NXDS ) 

' DO 1520 I = i ,NM9 
^ l T5 = NTH8B ( I ) 

1520 READ (5 , 1 603 ) < TH36 ( I »K ) , i<= 1 . NT3 ) 

DO 1 521 I = i. » NrVR 

NL=NLAve ( I ) 

1521 READ < 5 » 1 6^0 ) (LAM8 ( I »L ) . L= 1 i NL ) 

C THE INPUT FOR APPLYING THE PROGRAM TO A PARTICULAR VEHICLE BEGINS 
C WITH THE FOLLOWING READ STATEMENT AND CONTINUES THROUGH TEN 

c successive read statements 

9016 READ ( 5 » 9503 ) HEAD 

READ ( 5 * 9502 ) NCMPS , NMACHN , NOXS • NDX * NC 1 , NC2 * NC3 i NFCYCM , NTHFRR i NPUNCH 

RE AD ( 5 » 1600)THNR,TH3R,FAR,FNR 

READ ( 5 . 1600) XR4 , XR5.XR7.DX4 , DX5 . DX7 , XR l C 

READ ( 5 * 1600) (THFRP(N) . N=1 . NTHFRR ) 

READ ( 5 i 1600) (DORR ( J ) , J= 1 .NFCYCM ) 

READ (5 , 9502 ) (COMPNO ( I ) . I = 1 . NCMPS ) 

READ (5 i 9502 ) (NXS ( I ) . I =1 . NCMPS ) 

READ ( 5 . 1 600 ) ( DC (I ) « I = 1 . NCMPS ) 

READ (5, 1600 ) (MACHN (<) ,K=1 .NMACHN ) 

READ <5 » 1600 ) (XS <M ) , M= l , NOXS ) 

WR I TE(6t 9504 )HEAD 

WRITE (6 <90 04 ) NCMPS . NMACHN . NOXS »NDX *NC 1 » NC2 . NC3 , NFC YCM » NTHFRR » NPUNC 
1 H 

WR I TE ( 6 . 90 05 ) THNR , TH3R . FAR . FNR 

WRITE (6 .9307 )XR4. XR5, XR7.DX4 .DX5.DX7, XR1C 

WR I TE (6.9505 ) ( THFRR (N ) ,N= 1 ♦ NTHFRR ) 

WR I TE (6 . 9506 ) ( DDRR ( J ) » J = 1 . NFCYCM ) 

WRITE <6. 9010) CCCYPNO( I ) . DO ( I ) ,NXS( I ) , 1 = 1 .NCMPS) ' ' 

NM I = 0 
2000 LLL=0 
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17=1 
L 1 =0 
J= 1 
ND= 1 
L CN = 0 • 

FANS3=0 . 

MM = 0 
J5= 1 

THFR= THFRR ( j= ) 

DDR=DDRR ( J'5 ) 

NM I = NX 1 + 1 

IF ( NM I — NMACHN1 ) 90 1 5 « 90 1 5 »901 6 
9015 MR=MACHN (NM 1 ) 

IF(MR.GE • 15. ).VR=2 5. 

30 CO WRITE (6*901 1 )MR 
NBC1 =0 
KK2 = 0 

1000 LLL=LLL+ 1 

IF (LLL-NCMPS ) 536* 535 , 2000 
536 N = C0MPN0 (LLL ) 

DX = DO (LLL ) 

NX = NXS (LLL ) 

« 1 = K K 2 + 1 
KK2=KK1+NX-1 
MM = 0 

I 1 =NMACH (N ) 

GOT 0( 1 (2i3f4t5(6t1 ( 8) iN 
1 THMR1 =.C1 745*THNR 

LCN = DO (LLL )/(2.*(SIN(THNRl )/C0S ( THNRl ) ) ) 

DX=LCN 

1717 YAR=THNR 
K8=NXL 1 
J1 =1 

DO 1 1 1 3K= 1 ,NXL1 
1113 X (K ) =XL 1 (K ) 

CALL MACHXO ( I 1 i M * MR tLCNi XR 1 C » MM ,NDX« J 1 » DDR * I 7 ♦ THFR » NN 1 » 
XI «NX*ND«MWiv, j,h ,KK1 « KK2 » DX * LLL ) 

1 1 20 K3 = NTHN 1 ( I ) 

DO 1 1 1 9K = I ,K3 
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1119 YACK) =THN1 (I * K ) 

NL=NLAM 1(1) 

C01 1 12K=1 »NL 

1112 LAM (|< ) = LAM 1 ( 1 ,K ) 

GOTO 1 975 

2 LCN=FNR*DX 
DX=LCN 

917 Y AR = FNR 
<a=NXL2 
J1 = 1 

C042K= 1 * NXL2 
42 X ( K ) =XL2 ( K ) 

CALL MACHXDC I 1 , N * MR * LCN , XR 1 C * MM * NDX , J 1 * DDR , I 7 , THFR * NN 1 * 
XI * NX*ND*MMM, J,L1 i !<l< 1 ,KK2,DX,LLL ) 

47 :<3 = NFN2 < I ) 

NL = ML A M2 ( I ) 

D04 4l<= 1 » K3 

44 YA(!<) =FN2 ( I , K ) 

YAR=FNR 
D045K= 1 .ML 

45 LAM (K ) =LAM2 ( I ,K ) 

GOT 01975 

3 J 1 = 1 

CALL MACHXD ( I 1 , N , MR , LCN , XR 1 C , MM , NDX » J 1,. DDR *17. THFR * NN 1 * 
X I . NX . ND .VMM * J * L 1 *<K 1 . KK2.DX . LLL ) 

320 IF(NN1»EQ«1 ) GOT 03 39 
•Ml ( 1 ) = M A C H (3*1) 
vi (2 )=MAC.H (3*1+1 ) 

All (1 ) = A ( I ) 

A! 1 (2 )=A ( 1 + 1 ) 

CALL D I SCOT (MR, MR, Ml * A 1 1 , A 1 1 , - 1 0 * 2 * 0 * ANS ) 

340 CO 1 051 K = 1 , NX 

1051 AN3 (l< ) =8**3. 1 4*D0 (LLL ) *ANS* ( .5- XX (K ) )* ( 1 *-XX ( K ) )*XX (K ) 
NBC 1 = 1 

PANS3=AN3 (NX ) 

KK 3 = 0 

KK'<3 = !«2- 1 

CO 1 C52K =K!< 1 * KKK3 

KK3=KK3+ 1 
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WR I TE (6 . 9° 00 )XS (< ) , AN3 (KK3 ) 

IP (NPUNCH.EQ. 1 )G0TC555 
GOTO 1052 

555 PUNCH9000 , XS CK ) , AN3 ( i<i<3 ) 

1052 CONTINUE 
GOT 0 1 000 
330 ANS = A ( I ) 

GOT 034 0 

4 K9=NXD4 
J1 = 1 

61 D062K= 1 , NXD4 

62 X (K )=XD4 (K ) 

CALL MACHXD ( I 1 , N , MR , LCN , XR 1 C « MM * NDX . Jl . DDR iI7i THFR, NN 1 ♦ 

X I , NX,ND,M.vm, J,L 1 , << 1 i KK2 « DX , LLL ) 

I F ( NN 1 • EG • 1 ) G0T06 1 0 

1 F ( ( ( MACH (4,1 ) .C-E . 1 .35 ) .OR. (MACH (4 , 1 ) .LT. 1 .35 .AND. MACH (4 ,2 ) .GT. 
XI .35 ) ) .AND. (MR.LE. 1 .35 ) JG0T063 
610 !<3 = NFN4 ( I ) 

NL =NL A M4 ( I ) 

C064K = 1 . K3 

64 YA (K )=FN4 ( I , '< ) 

YAR=FNR 
D065K= 1 » NL 

65 LAM (K ) =LAM4 ( I , K ) 

GOT O 1 075 

63 WRITE (6. 9003 ) 

GOT 01000 

5 YAR=THNR 
K 8 = NX D 5 
Jl = 1 

C052K = 1 . NXD5 

52 X(K ) = XD5 (K ) 

CALL MACHXD ( I 1 , N , MR » LCN . XR 1 C , MM . NDX . J 1 .DDR, 17, THFR , NN 1 , . .. 

XI ,NX,ND,MMM, J, LI 1 ,KK2,DX, LLL ) 

520 K2=NTHN5 ( I ) 

NL=NLAM5 ( I ) 

D05532K = 1 ,'<3 
5532 YA (K ) =THN5 < I ,< ) 
w D05533K=1,NL 
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tn 

•*“ 5533 LAM (K ) =LAM5 ( I ,K ) 

1075 CALL INTP(N.NX,YAR,NL,K8,NN1 , J1 , I ♦ LCN , FNR , F ANS2 ♦ FANS3 » MR , NC 1 ♦ 

XX R4 , DX4 .NDX.NC2, XR5 »DX5i K< 1 , KK2 ,LLL , MM , NBC 1 .NPUNCH) 
IF(NN1.E0.1.0R.J1*GT.2 )G0T0775 
1013 IF (N.EQ. 1 )G0T01 120 
IF(N.EQ.2)G0T047 

IF (N. EQ. 5. OR. N. EG. 6. OR. N.EQ. 7 JGOT0520 
I F (N.EQ. 4 )G0T061 0 
IF (N.EG.9)G0T01 304 
775 IF (N.EQ. 6 )G0T08276 
IF (N.EQ. 7 JGOTOI 750 
GOT 0 1 000 

1750 MR = MR7 
NN1=NN707 
1 = 17 

I F ( ( ( DDR-DDR7 (1,1 ) ) . LT • - • 1 E-7 ) • AND .L 1 . EQ • 1 ) GOTO 1 059 
IF ( ( (DDR-DDR7 ( I . 1 ) ) . LT . - • 1 E-7 ) .AND. LI . EQ • 2 • AND • ( (THFR-THF7 ( I » J , 1 
X) ) • LT • — • 1 E— 7 ) ) GOTO 1751 
I F ( (THFR-THF7 ( I ♦ J , 1 ) ) . L T. - . 1 E-7 ) GOTO 1 751 
I F ( (DDR-DDP7 (I , NO ) ) . GT . . 1 E-7 ) GOTO 1 060 

1751 N3C= 1 

DO 1 061 K2 = 1 .NX 
XX ( <2 ) =XX7 (K2 ) 

1061 A A 1 ( K2 ) =A4 (l<2 ) 

THFR=THF7 ( I , J, 1 ) 

GOTOl 062 

1059 THNR=THNR 1 

DO 1 067K2= 1 ,NX 
1067 A 1 ( 1 ,K2 )=A4 C<2 ) 

GOT 0 1 053 

1060 DO 1 065K2 = 1 ,NX 

A 1 (2.K2 )=A4(K2 ) 

1065 XX (K2 ) =XX7 (K2 ) 

GOTO 1 053 

8276 D0772SK2= 1 ,NX 
XX ( K2 ) =XX6 (K2 ) 

7728 A 1 ( 1 ,K2 ) =A4 (<2 ) 

NBC= 1 
G0T08275 


APPENDIX A 



6 MM = 6 

THFR1 = % 0 1 745*THFR 
NN2 = 0 
Jl = l 

I F ( (FAR-4. ).GE#-# 1E-7)FAR=4. 

CALL MACHXD ( I 1 i N , i*R *LCN * XR 1 C * MM *NDX* J 1 * DDR *17* THFR * NN 1 * 
XI * NX * NO * MMN- < J,L 1 *<<1 *KK2*DX,LLL ) 

81 03 LI =1 
L2= 1 
L 3= 1 
L4= 1 
E 1 =0 • 

E2 = 0 • 

E3 = 0 • 

E4 = 0 • 

NN3 = 0 
NN4 =0 
NBC = 0 

250 NTN=NTHN6 ( I ) 

D0251 J= 1 * NTN 

I F ( ABS ( THNR-THN6 ( I * J ) ) . LE • • 1 E-7 )G0T0252 

251 CONTINUE 

IF (THNR*LT.THN6 (1*1) )G0T0253 
D091 35J=2*NTN 

IF(THNR.LT.THN6 ( I * J ) )G0T0254 
9135 CONTINUE 
J=NTN- 1 
GOT 08254 

253 J=1 

GOT 08254 

254 J= J-l 
G0T08254 

252 J=J 
NN2 = 1 

8254 NTF = NTHF6 ( I * J ) 

D08257K = 1 * NTF 

I F (ABS ( THFR-THF6 ( I * J * K ) ) . LE . • 1 E-7 ) G0T0825Q 
8257 CONTINUE 

8256 I F ( THFR.LT.THF6 ( I * J * 1 ) ) GOT 08260 


ui 
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D0826QK=2<NTF 

I F ( THFR.LT .THF6 (I ) GOT 0 8255 

8268 CONTINUE 

k=ntf-i 
GOT 0826 1 
8255 K=K- 1 

GOT 0826 \ 

8258 ]< = K 

NN3 = I 
Got 0826 1 

8260 D08270l<2= 1 ,NX 
8270 XX6 (K2 ) =XX ( K2 ) 

8273 NN606 = NN 1 

THNRS=THNR 
THNR = THN6 ( I f J ) 

THFRS=TH- R 
16=1 
MR6=MR 

WR=MACH (6^1) 

I 1 =NMACH (5 ) 

N = 5 
GOT 05 
8275 1=16 

NN1 =NN606 

THNP = THNRS 

T HFR = THF6 (UJ,1 ) 

MR=MR6 
LI =2 
K=1 

8261 IF(< # F0.1 ) G0TQ8755 
LL1 = l 

LL 2 = O 
LL 3 = 1 
LL 4 = 0 
!< K = K — 1 

0 0 1 999L = I U<K 

LL 1 =LL1 +NLAM6 (hJ»L) 

1999 LL3=LL3+NFA6 ( I . Ji L ) 

002002L = 1 •< 
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LL2=LL2+NLAw6 ( I iJiL ) 

2002 LL4=LL4 + NFA6 ( I , J,L) 

G0T080 02 
8755 LL 1=1 

LL2 = NLA^6 < I , J, 1 ) 

LL3= 1 

LL4 =N P A 6 ( I ( J * 1 > 

80 02 NL =0 

D08003L=LL1 .LL2 
\'L=NL+1 

8003 LAM(NL)=LAV6( I ,J«L) 

NF = C 

D08004L=LL3,LL4 
NF=N r + 1 

FA ( M p ) =FA6 ( I i JiL ) 

8004 K 1 1 (NF)=K1 ( I « J « |_ ) 

C08C 1 K2 = 1 » NX 

XX 1 =XX ( <2 ) 

CALL D I SCOT (FAR, XXI (FA ( LAM « XD6 » 1 1 , NL , NXD6 , A NS 1 ) 
801 A 1 ( L 1 ,K2 )=ANS1 
D0890 1 L= 1 (NF 

IF (ABS (FAR-FA (L ) ) .LE. . 1 E-7 )GOT08903 

8901 CONTINUE 

I F (FAR.LT.FA ( l ) ) GOTO8902 
C08 1 L =2 * NF 

IF (FAR.LT.FA (L ) )G0T089Q4 
81 CONTINUE 
L=NF— 1 
GOT08905 

8902 L=1 

GOT 08905 

8904 L = L- 1 
GOT 08905 

8903 L=L 
NN4 = 1 

8905 IF (NN4.EQ. 1 )G0T08906 

I p (K1 1 (L ) .EQ. 1 .A ND.X1 1 (L+l ) .FQ. 1 ) GOT 0259 
IF ( K 1 1 (L ) .EQ.O • A NO • K 1 1 (L+l ) .EO. 0 ) GOT 081 5 
I F (K 1 1 ( L ) • EO • O • AND « K 1 1 (L+l ) .EO. 1 ) GOT 08909 


cn 
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I F ( K 1 1 (L).EQ.l .AND. K1 1 (L+l ) .EQ.O ) GOT 089 10 
8906 IF (K1 1 (L) «EQ.0)G0T08l5 
GOT 0259 

8909 K 6=L 

802 K5=NL 
K3 = NF 
G0T081 4 

8910 K6=L+1 
GOT0802 

81^ El =F 1 +1 . 

MN 1 1 =L1 
GOT 0259 

814 D08 1 6NN = K6 « <5 » K3 
816 LAM (NN ) = 1 . 1 *LAM (NN ) 

XXI =XX ( 1 ) 

CALL D I SCOT (FAR. XXI . F A . LA.M ♦ XD6 . 1 1 .ML .NXD6.ANS2 ) 
I F ( ABS ( ( A 1 (LI . 1 ) — ANS2 ) / A 1 ( L 1 . 1 ) ) . LE . . 02 1G0T0259 
FI =E1 +1 . 

MN1 1 =L1 

259 I F (NN3.EQ. 1 1GOT0260 
LI =L1 +1 

I F (LI -2 ) 2 583 ,2^88.261 
2588 K=K+1 

G0T0826 1 

261 IF (NBC.EQ. 1 )G0T08Q7 
THF ( 1 ) = THF 6 ( I . J , K— 1 ) 

THF (2 )=THF6( I . J.K ) 

GOT 0808 

807 THF (1 ) =0 % 

THF (2 ) = THF 6 ( I . J , 1 ) 

THFR=THFRS 

808 C0803K2 = 1 .NX 
All (1 )=A1 (1, <2 ) 

All ( 2 ) = A 1 ( 2 ,K2 ) 

CALL 01 SCOT ( THFQ, THFR, T HF . A 1 1 , A 1 1 . - 1 0 . 3 . 0 . ANS 1 ) 

803 A2 (L2.K2 ) =ANS1 

IF (El »LE. • 1 E— 7 ) C0TO262 
I F ( (El-1 . ) .LE. . IE-7 1GOTO805 
8C6 E2=E2+1 • 
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MM1 1 = L2 
G0T0262 

805 A1 1 (MNii 1 j = l . ]»A1 (MN1 1 ( 1 ) 

CALL D I SCOT ( TH^R, THFR i T HF i A 1 1 , A 1 1 i-10i2i0.ANS2) 

I - ( A5S ( ( A2 (L2 , 1 ) - AN S2 )/A2 <L2 • 1 ) ) .LE. .02 JG0T0262 
GOT 0806 

260 C0263<2= 1 ,NX 

263 A2-(L2.K2 ) = Al (LI »<2) 

E2 = E1 

262 I F ( NN2 • EO • ! ) GOT 02 65 
L2 = L2 + 1 

IF (L?-2 ) 266 *266 » 2 67 

266 Ll=l 

J = J+ 1 • & 

FI =0. 

GOT 08254 

267 TN ( 1 ) =THN6 ( I »J-1 ) 

TN (2 ) = THN6 ( I i J ) 

D08077K2= 1 (NX 
A22 (1 )=A2 ( 1 ,<2 ) 

A22 (2 )=A2 (2 (K2 ) 

CALL D I SCOT ( THNR * THNR i TN i A22 » A22 i-10(2*0i ANS 1 ) 
8077 A3 (L3»K2 }=ANS1 

IF (E2.LE. .1 E-7 ) GOT0268 
I F ( ( E2- 1 . ) .LE. . 1 E-7 )GOT08808 
809 E3=E3+1 . 

MN1 1 =L3 
G0T0268 

8808 A22(MN1 1 ) = 1.1*A2(VN1 1 »1 ) 

CALL D I SCOT ( THNR , THNR , TN , A2 2 » A22 (-10(2«0»ANS2) 
IF ( ABS ( (A3 ( L3 ( 1 ) -ANS2 )/A3 <L3 i 1 ) ) . LE . . 02 ) G0T0268 
G0T0809 

265' D027n<2=l , NX 

270 A3 (L3,K2 )=A2 (L2 ,K2 ) 

E3 = E2 

268 IF(NN1 .EO.l 1G0T0271 
L3=L3+ 1 

IF (L3-2)2 72 »272»273 
272 1=1+1 


I 
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L 1 = 1 

L2=1 flu v:3~c 

E 1 = 0 • /U/V'S-C' 

c 2 = 0. 

G0T0250 

271 D0274!<2= 1 ,NX 

274 A4 (K2 )=A3 (L3 «K2 ) 

E4 = E3 
GOT 0277 

273 Ml ( 1 )=MACH ( 6 * I - 1 ) 

Ml (2 ) = M A C H ( 6 i I ) 

D086‘72l<2= 1 .NX 
A33 ( 1 ) =A3 ( 1 ,K2 ) 

A33 (2 ) =A3 (2 , K2 ) 

CALL DI SCOT (MR ,MR,M 1 , A33 , A33 * - 1 0 i 2 » 0 * ANS 1 ) 

8672 A4(«2)=ANS1 

I f (E3.LE . . 1 E-7 ) G0T0277 

IF( ( E3- 1 . ) .LE. . 1 E-7 ) G0T0278 

E4 = E4+ 1 . 

G0T0277 

278 A33 ( MN 1 1 ) = 1 . 1 *A3 ( MN 1 1 » 1 ) 

CALL D I SCOT ( v.R * MR , M 1 , A3 3 * A3 3 , - 1 0 * 2 » 0 » ANS2 ) 

I F ( ABS ( ( A4 ( 1 ) -ANS2 ) /A4 ( 1 ) ) • LE . • 02 ) GOT02 77 
E4=E4+ 1 . 

277 K K 3 = 0 

D08277K2=KK1 %KK2 
KI<3=KK3+ 1 

CX = DO (LLL )*< 1 .+2.* XX (Kl<3 )*S I N ( THFR1 ) /COS (THFR1 ) ) 
ANS = DX*A4 ( '<K3 ) 

WRITE (6 *9000 )XS (K2 ) , ANS 
IF (NPUNCH.EQ. 1 ) G0T0556 
G0T08277 

556 PUNCH9000 * XS (K2 ) * ANS 
8277 CONTINUE 
GOTO 1 000 
7 M M = 7 
1049 NBC = 0 
LI =1 
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L2=l 
\!\i2 = ^ 

CALL MACHXD ( II «N« : ^R * LCN * XR 1 Ci^-v *NDX, J1 * DDR , I 7 * THFR , NN 1 , 
XI iNX,ND» >1MV , J ,L 1 * 1 * Kl< 2i DX , LLL ) 

1 021 ND=NDDR7 ( I ) 

IF( (DDR-DDR7 (I*1)).LT.-. 1 E-7 ) GOTO 1 0 57 
D01077J=1 «ND 

I F ( ABS ( DDR-DDR7 ( I « J ) ) • LE • • 1 E-7 ) GOTO 1 02 7 
1077 CONTINUE 

I F ( ND • GE • 2 ) GOT 01776 
J= 1 

GO F01035 

1776 DO 1 02 8 J = 2 * ND 

IF (DDR.LT.DDR7 ( N J ; ) GOTO 1011 

1028 CONTINUE 

1029 J = NR 
GOTOl 035 

1011 J=J-1 
GOTOl 035 

1027 J = J 

NN2 = 1 

1035 I F ( ( THFR- THF7 ( I iJil ) ) • L T • - • 1 E-7 ) GOTO 1 050 
GOT 01062 
1050 THFR 1 =TH r P 
I 7= I 

DO P 1 2K = 1 4 NX 

1012 XX7 ( K ) = XX (!< ) 

1063 NN707=NN 1 

MR7=MR 

MR=MACH (7,1 ) 

I 1 =NMACH ( 5 ) 

N = 5 
GOT 05 

1057 THNR 1 =THNR 
THNR = THFR 

I 7= I 

1058 DO 1 775K = 1 ,NX 
1775 XX7 ( K ) = XX (K ) 


a 


APPENDIX A 



GOTO 1 063 

1 062 NTF = NTHF7 < I , J) 

NL = NLAM7 ( I « J ) 

DO 1 064L = 1 .NTT 
1064 THF (L ) =THF7 ( I , J ,L ) 

DO 1 079|_= 1 iNL 
1079 LAM (L ) =LAM7 ( I , J ,L ) 

D01055I<2=1 ♦NX 
XX1=XX(I<2) 

CALL D I SCOT (THFR. XXI , THF » LAM , XD7 . 1 1 , NL i NXD7 , ANS ) 

1 055 A1 (LI .K2 ) =ANS 

I F ( NBC • EQ • 1 ) GOTO 1089 
GOTO 1 053 
1089 THFR=THFR 1 
THF ( 1 ) =0 . 

THF ( 2 ) = THF7 ( I , J,1 ) 

D0777!<2=1 .NX 
AA2 ( 1 ) = AA 1 (!<2 ) 

A A2 (2 ) = A 1 (LI . K2 ) 

CALL D I SCOT (' THFR . THFR .THF.AA2.AA2.-10.2.0.ANS) 

777 A 1 (LI .K2 ) =ANS 
1053 IF (NN2.EQ. 1 1G0T01 086 
LI =L1 +1 

IF (LI -2 ) 1 045 . 1 045. 1 088 

1 045 I r ( ( DDR-DDR" 7 ( I , 1 ) ) . L T - - • 1 5-7 ) GOTO 1 046 

I F ( ( ( DDR-D0R7 ( I ,ND) V.GT * • 1 E-7 ) .AND. ( ( f H F R - T HF 7" ( I , J ,1 ) V'.LT. 
X— • 1 E-7 ) ) GOT 01046 

I F ( (DDR-DDR7 ( I . NO ) ) .GT . . 1 E-7 ) GOTO 1 050 ’ 

J = J+1 
GOTO 1035 
1048 J=1 

GOTO 10 35 

1 046 DO 1 0 4 7 1< 2 = 1 .NX 
1047 A 1 (LI . K2 >=AA1 (K 2 ) 

1 088 I F ( (DDR-DDR 7 (I , l ) ) . LT . - . 1 E-7 ) GOTO 1 068 
I F ( (DDR-D.DR7 ( I . NO ) ) .GT. . 1 E-7 ) GOTO 10 70 
D2 ( 1 ) =DDR7 ( I , J- 1 ) 

D2 (2 ) =DDR7 ( I , J ) 
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1069 D01C9S'<2=1 .NX 

All (1 ) =A1 ( 1 ,K2 ) 

All (2 ) = A 1 (2 , «2 ) 

CALL 0 I SCOT ( DDR .DD5 . D2 .A1 1 , A 1 1 .-10.2.0. ANS ) 

1098 A2 (L2.<2 ) =ANS 
G0T01 085 

1068 D2 ( 1 >=0. 

D2 (2 ) = DDR7 (1,1) 

GOTO 1 069 

107^ D2 ( 1 ) =DDR7 (I ,NO ) 

D2 (2 ) = 1 . 

G0T01 069 

1 086 COl 087K2= 1 ,NX 
1087 A2(L2,K2)=Al (1 ,K2) 

1085 IF(NM1»EG«1 ) GOTO 1091 
L2=L2+1 

I F (L2-2 ) 1 092 , 1 092 , 1 093 

1092 LI =1 
NN2 = 0 
1 = 1 + 1 
G0T01 02 1 

1093 Ml ( 1 )=VACH<7, 1-1 ) 

Ml (2)=MACH(7,I ) 

DO 1 0 9 4 K 2 = 1 ,NX 

A 1 2 ( 1 ) =A2 ( 1 ,K2 ) 

A 1 2 (2 )=A2 ( 2 , K2 ) 

CALL D I SCOT ( MR ,MR ,M 1 , A1 2, A1 2 ,-l 0 .2 , 0, ANS ) 

1094 A4(K2)=ANS 
G0T01 097 

1091 KK3=0 

DO 1 09 <?i< 2 = I<K 1 ,K!<2 
KK3=K<3+1 

1099 A 4 (KK3 ) = A 2 ( 1 ,KK3) 

1097 KK3=0 

DO 1 096K2 = K'< 1 ,I<K2 
<K3=KK3+1 

IF ( J5.EQ.NFCYCM ) GOTO 1250 
GOTOl 231 

1250 IF (NC3.E0. 1 JG0T0122B 


i 
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1231 ANS=A4 (KK3)*DX 
GOTO 1 227 

1228 IF (ABS <XS (K2 1-XR7 )-.LE'» . 1 E-7 IGOTOl 229 
GOTO 1232 

1230 ANS=A4 (KK3 )*DX7 
1227 WR I TE ( 6 » 9000 )XS (K2 ) iANS 
IF (NPUNCH.EQ.l 1G0T0557 
GOTO 1 096 

557 PUNCH9C00.XS <K2 ) ♦ ANS 
GOTO 1 096 

1229 ANS = A4 (l<l<3 )*DX 

WRITE (6.9000 ) X S ( K 2 ) tflNS 

IF (NPUNCH.EQ. 1 )G0T055S 
G0T01 230 

558 PUNCH9000 , XS (K2 ) ♦ ANS 
GOT 0 1 230 

1232 IF C (NDX.EO.O.AMD.XS (K2 ) .GT.XR7) .OR. (NDX.EQ.l . AND.XS(K2) .LT.XR7) ) 
XGOTOl 230 

GOTO 1 2 3 1 
1096 CONT I NUF 

IF (NFCYCM.GT.l )G0T01292 
GOTO 1 000 
1292 J5=JF+i 

IF ( J5.GT. NFC VC M) GOTO 10 00 
THFR=THFRR ( J5 ) 

DDR = DDRR ( J5 ) 

F AR = 4 . 5 
G0T01 000 
8 YAR=THBR 
'< 8 = N X D 8 
J1 = 1 

DO 1 30 1 K = 1 ,NXD8 
1301 X(K)=XD8(K) 

CALL MACHXD ( I 1 . N. MR . LCN « XR 1 C . MM , NDX . J 1 % DDR iI7i THFR «NN1 i 
X I * NX.ND.MMM , J.Ll »«1 . KK2 0XiLLL ) 

I F ( NN 1 • EQ • 1 ) GOT 0 1 304 

I F ( ( (MACH (8 , 1 ) . GE. 1 . 35 ) .OR . (MACH (8 . 1 ) .LT . 1 • 35 . AND. MACH ( 8 ,2 ) .GT . 
XI .35 ) ) .AND. (MR.LE. 1 .35 ) ) GOT 09 1 
1 304 !<3 = NTHB8 ( I ) 
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DO 1 '302K = 1 .<3 
1302 Y A ( K ) = T HB8 (I » K ) 

ML=NLAM8 ( I ) 

DOl 303K= 1 ,NL 
13 03 LA iv (K ) = LAM8 ( I , K ) 

GOTO 1 075 

91 V/P I TEC 6. 90 06 ) 

GOTO 1 000 

1600 FORMAT ( 1 CF7, 3 ) 

9000 FORMAT ( 2E 1 6 • 8 ) 

9001 FORMAT ( 1 X66HTHE DLAM5DA VALUES r 0R THE FRUSTUM ARE ESTIMATED FOR T 
1 H I S MACH NO. ) 

9005 FORMAT C 1 X5HTHNR = FQ. 3, 3X5HTH6R = F8 . 3 , 2X4 HFAR = F8 . 3 , 4 X4HFNR= F8 . 3 ) 

90 07 FORMAT ( 1 X4 HXR4 = F9 . 3 OX 4HXR5= F9 . 3 . 2X4 HXR7 = F8 . 3 , 4X4HDX4 = F8 . 3 , 

X3X4HPX5 = F8.3 »3XAHDX7 = F7 .3* 4X5HXR1 C = F3.3 ) 

9502 FORMAT (1415) 

9503 FORMAT ( 1X1 1 A6 ) 

9504 FORMATC 1H1 , 1 1A6/40X16HINPUT PARAMETERS) 

9004 FORMAT ( 1 X6HNCMPS= I 5 . 5X7HNMACHN= 15. 3X5HN0XS= I 5 . 6X4HNDX= I 5 . 6X4HNC 1 = I 
1 2 . 3X4HNC2= I 2 ♦ 3X4HNC3= I 2 ,3X7HNFCYCM= I 2 , 3X7HNTHFRR= 1 2 , 3X7HNPUNCH= I 2 ) 
9010 FORMAT ( 1X/1X9HC0MPN0 ( I )5X5HD0 ( I )5X6HNXS ( I ) / ( I 5 . 7XF8 .3.4X15 ) ) 

901 1 FORMAT ( 1 X/ 1 X5HMACH = F7 . 3//6X7HST AT I 0N9X7HDLA MB DA ) 

9003 FORMATC 1X101HNO FIXED INPUT DATA FOR CYLINDER FOLLOWING OGIVE IN S 
1UBSONIC AND TRANSONIC RANGE ARE INPUT TO PROGRAM) 

9006 FORMAT ( 1 XI 09HNO FIXED INPUT DATA FOR BOATTAIL FOLLOWING LONG CYLIN 

1 DER IN SUBSONIC AND TRANSONIC RANGE ARE INPUT TO PROGRAM) 

9505 FORMAT ( 1X/5X8HTHFRR ( I )/( 1 X5FS.3 ) ) 

9506 FORMATC 1X/5X7HD0RRC I )/( 1X5FQ.3 ) ) 

END 


Oi 

tn 
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SUBROUTINE MACHXD ( I 1 , N * MR * LCN * XR 1 C * MM , NDX , J I , DDR , I 7 , 

XTHFRiNNl , I J,L1 ,<.<1 ,KK2*DXiLLL) 

C CALCULATION of MACH NUMBERS USED IN INTERPOLATION AND CONVERSION OF 
C STATIONS INTO X/DO OR X/L VALUES 

REAL MACH, Ml ,LCN , MACHN , LAM 1 , L.AM2 , L AM5 , L AM6 , L AM7 , LAMQ , lAM , MR 
COMMON NMACH (8 ) ,NTHN1 ( 1 A ) , NLAM1 ( 1 4 ) , MACH ( 8 , 1 6 ) * THN 1 ( 1 4 , 1 7 ) , XL1 (6 ) , 
1 LAM1 ( 1 4 i 1 02 ) ,NLAM2 ( 7 ) , NFN2 ( 7 ) , XL2 ( 2 1 ) , FN2 ( 7, 5 ) ,LAM2 ( 7i 1 05 ) , A ( 12 ) , 
2NLAM4 ( 5 ) , NFN4 ( 5 ) * XD4 ( 14 ) , FN4 ( 5 , 3 ) , LAM4 ( 5 , 42 ) , NLAM5 (16) , N THN5 (16) , 
3XD5 (24 ) , THN5 ( 1 6 ,9 ) , LAMB ( 1 5, 2 1 6 ) ,NTHN6 (3 ) ,NTHF6 (3 ,3 ) , 

4NFA6 < 8, 3 , 4 ) ,NLAM6 ( 8 , 3 , 4 ) , NT0TL6 <8 , 3 ) ,NTOTFA (8 ,3 ) ,THN6 (8,3) , 

5XD6 (18) , THF6 (8,3,4)*LAM6(8,3,288),FA6(8,3,16) ,I<1 (8,3,16), 

6NDDR7 (8) ,NTHF7(8,3) ,NLAM7(8,3) , DDR7 ( 8 , 3 ) , ThF 7 ( 6 , 3 , 4 ) , L.A M 7 ( 8 , 3 , 96 ) , 
7XD7 (24 ) , NT HR 8 (8),THb8(S,5) *XD8( 12) , LA.M8 (8 ,60 ) , NLAM8 ( 8 ) 

COMMON HEAD ( 1 1 ) ,XX( 100 ) ,X (25 ) , LAM (225) ,NXS( 14 ) ,XS(350 ) *D0( 14) , 

1 A 1 ( 2 , 1 00 ) , YA ( 1 7 ) , TN ( 2 ) , MACHN (15) 

COMMON A2 ( 2 , 100 ) , A3 (2 , 1 00 ) , M 1 (2), All (2),A4(100),AA1 (100), 

1 AN3 (5 ) , XX 6(100) ,FA ( 4 ) , <1 1 (4) , THF ( 4 ) , A22 < 2 ) ,A33(2 ) ,XX7( 100) ,AA2(2) , 
2D2(2) , A 1 2 ( 2 ) ,COMPNO ( 1 4 ) , A44 ( 2 , 1 00 ) , THFRR ( 5 ) , DDRR (5 ) 

IOC NN 1 = 0 

DO 1 0 1 I = 1 , I 1 

IF ( ABS ( MACH (N, I ) -MR ) . LE • • 1 E-7 ) GOTO 1 02 

101 CONTINUE 
GOTO 1 03 

102 NN 1 = 1 
I = I 

G0T01 04 

103 I F (MR.LT.MACH (N , 1 ) ) GOTO 1 05 
IF (MR.GT.MACn (N, I 1 ) ) GOTO 106 
CO 1 07 1 = 1,11 

I F ( MR-MACH ( N , I ) >108,108,107 

107 CONTINUE 

108 1=1-1 
GOTO 104 

105 1=1 
GOTO 1 04 

106 1=11—1 

1^4 IP - (N.E0.5 )GOTQ41 1 3 
G0T041 1 2 

4113 IP (MM.EQ.6 )G0T041 15 
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'I F (MM, eg. 7. A\JD. ( ( D0R-DDR7 ( 17,1 ) ) . L T . - . 1 E-7 > . AND . (LI »£Q» 1 ) 

X ) GOT 04 1 2 1 

IF (MM. EQ .7* AND. ( < DDR-DDR7 ( 17,1 ) ) . L T . - . 1 E-7 ) . A ND . L 1 .E 0.2. AND. 
X( ( THFR- THF7 ( I 7 , J , 1 > ) . LT . - . 1 E - 7 ) ) G0T04 1 2 2 
I F < (MM.EQ.7) .AND. ( ( THFR-THF7 ( 17, J, 1 ) ).LT.-.lE-7) 1G0T04122 
I F ( (MM.E0.7 ) .AND. ( ( DDR -DDR 7 ( I 7 , ND ) ) . GT . . 1 E-7 ) )G0T041 22 " 

GOT C A 1 1 2 
4115 N = 6 

XR 1 = XR 1 C . 

GOT04109 

4121 N = 7 
G0T041 12 

4122 N = 7 

XR 1 =XR 1 C 
G0T041 09 

4112 IF (N.EQ. 1 .0R.N.EQ.2 )G0T04104 
XR 1 =XS ( KK 1 ) 

4109 IF(NDX.EQ.l )GOT04107 
KK3 = 0 

DQ4 1 06K=K!<1 , <K2 
KK3=K<3+1 

4106 XX(K<3)= (XS(K)-XRl )/DX 
GOTO 110 

4107 KK3 = 0 

D04 1 08K = K< 1 , K!<2 
KK3=KK3+ 1 

4108 XX (KK3 ) = ( XR1-X5 (K ) )/DX 
GOTO 110 

4104 IF (NDX.EQ.O ) XR 1 =XS ( KK 2 ) -LCN 
IFtNDX.EQ.l )XR1 =XS(K!<2 )+LCN 
GOT04109 

110 RETURN 

FND ' 
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03 

CO 


SUBROUTINE INTP ( N , NX , YAR , NL , K8‘, NNf, 0 Y, I , LCN <FNR<FANS2< FA NS 3 » MR « 

XNC 1 , XR4 i 0X4 » NDX , NC2 » XR5 , 0X5 , KK 1 ,KK2 *LLL , MmV'NBC Y, N PUNCH j 

C INTERPOLATION ROUTINE FOR - COMPONENTS “1 ♦" 2 V 4 T 5 * AND '§ ' 

REAL' WCHTmR ,'mTYLCN 7i^'CffNTL AM 1 YlA M2TE AM 4TlA M 5Tl A M6 , LA M Y» L AM 87lAM 
COMMON NMACH (8 ) , NTHN j. ( 14), NLAMl'Vi 4 )TmACH ('8Yr6~)TT HNY( 14 /1 7 )~, X'LTf6 
T'LAMTc 14, 102 f, NLA M2 (7'T7 nfn2TT >1 XL272TT,' FN 2 1 77 5~) 7U A M 8 T T % T05 F, A ( 1 2T,“ 
' 2NLAM4 (5 ) YNFN4 ( 5 ) , XD4 (Y'4 ) '» FN4 (5,3) , LAmATY, 42T ,"NL A M 5 (Y6 F, nTHN" 5 T 16T, _ 

3X05(24 ) , THN5 (16,9) , L AM 5 ( 1 6 » 2 1 6 ) ,NTHN6~( 8 ) , NTHF6 (8,3) , 

' 4NFA6 ( 8 , 3% 4 > , NL AM6 ( 8 ,3,4), NT 0 TL6 ( 8 , 3") » NTOTFAT 8 T 3 j , _ THN6 ( 8, 3 ) Y 
YYD6 (18 ) , THF6 (8 ,3,4 ) , LAY 16 T 8 ,3,288) ,FA6 ('§73 , 1 6 )TK1 (8, 3~,”l6T, 

6NDDR7 ( 8 ) , NTHF7 (8,3), NLA M7 ( 8", 3 ) , DDR7 7873 ~j YtHF 7(8,3 , 4 ) V LXmT( 8, 3 , 96T7 
7X07“ ( 2A~) 7 N T H B8 T 8 ) ,YHB8 ( 8 TYTA XD8 ( 1 2 )YlAm8‘( 8,60) , NLAM8 ( 8 ) 

COMMON 'iHEADY' ri )“, XX flOO ) , X (25“) , LAM (“225 ) , NXS ( 14 ) ,X~S (350 ) , DO ( 1 4 ) , 

1 A1 (2, 1 00 ) , YA ( 1 7 ) , TN (2 ) , MACHN ( 1 5 ) ' 

COMM ON A2 ( 2 , 1 0 0 ) , A3 ( 2 , 10 Ol , MT ( 2 ) , A 1 1 ( 2l 7A4 '(To 0 ) TA AT ATO Ol“, 

T AN3(5Y VXX6 (lOO ) ,FA ( 4 ), K 1 1 (4 ) A THFT 4 1/A221 217 A 3 3 (HTf, X X 7 (_1 0(7 ) , A A2 ( 2) , 

202(2), A 12 (2 ) ,COI^ PNO ( i'4 ) , A44 ( 2,"l 00) , f HFRR (' 5 ' T/DDRR (5 j" 

001 6 76l< 2 = 1 ,'NX 

XX I =XX ( K2 j " 

CALL OISCOT ( YAR,XX1 , Y A , LAM 7 X , 1 1 , NL , K8 , ANS ) 

1076 A44 < J1 ,K2 ) =ANS 

1 0 1 1 _ I FTnnF/E^Q • ll GOTOl 015 

7)1 =J1 +1 

pr ( j j - 2 J 1 cT3 , f 6 1 3 , TO T4 " _ 

T 0 T 3 i"=T+T ' 

RETURN 

T0“14 If (MM.'EQ.6,0R7^M.EQ.'7”j"N=‘5' — ' ' 

M 1 < i ) =MACH (N , I - 1 ) _ 

Ml (2 ) = MACH (N , I ) “ 

'0051 08K2= 1 ,NX " 

Al l (1 ) = A44 ( 1 ,K2 ) 

All (2)=A44(2,K2) " ' “ 

CALL D I SCOT ( MR , MR , M 1, A'ff , A 1 1 ,-10,2,0, ANS) “ 

'5108 A4 ( K2 ) = ANS ' — 

GOT05109 ' . .. - 

" i01'5 C'0'101 6K2 = 1 .NX ^ " ' 

10 1 6" A‘4TK2 j=A"44 ( 1 ,K2 ) '' 

51 09 'kl<'3 = 0 ' * " ' ' 
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'IF(MM.EQ.6»0R.MM.EQ.7)N = MM 

D05i 1 0K2=KK 1,K'<2 

KK3=KK3+ 1 
I f(n.EG 7 1 ) G0TC9 i 6 

IrlN.EQ _ .2 JGOTOI'OI 7 " ' ' " 

IF (N.EQ.4.0R.N.EQ.5 JG0T01020 
I F ( N • EQ • 6 • OR • N • EQ • 7 ) GOT 05 110 
I F ( N • EQ • 8 ) GOTO 1 305 
916 THNR1 = • 0 1 745*YAR 

DX=2V*(S IN (“THNR1 ) /COS ( THNRl) ) *XX (KK3 ) *LCN 

G 0 T 01 q 19 

1017 CX = 2 • *D0 (LLL )* ( .25-FNR*'*2+SGRT ( ( FNR**2 + . 25 ) *.*2- ( XX ( <l<3 ) *FNR-FNR 

” T )**2) > 

GOTO 10 19 

1305“ fHBRl = .01745*YAR 

"D X = D 0 (LLL )* ( 1 . -2.*XX ( KK3 )*S IN (TH6R1 ) /COS CThBRl ) ) 

GOTO 1019 

T620 DX=D0 (LLL ) 

"10T9” ANS = DX*A4 "(KK3 ) 

rFTN.EO. 1 1G0T0774' 

IF ( N • EQ • 2 ) G0T0776 

f'F'('N.EQV4 _ )G'0T0779 

TF(NVE6V5 ") GOTO 778 

I F ( N»EG • 8 ) G0TC762 ' ~ ' ~ 

774 I F ( NBC 1 . E G . 1 . A ND . KK3 . EQ . 1 ) ANS= ( ANS+FANS3 ) /2 “ 

WRITE ( 6 « 90 00 ) X STI< 2 ) VANS 

I F ( NPUNCH , EQ . 1 ) G 0 T 0 5 5 5 " ' 

G0T05 f T6 " ~~ 

~ *553“ PURc oo 6^ 5< 57K2 T >' Tans " ‘ 

'GOT 05 ITO 

“776"TF ( NBCT YEQ . l '.AND'. KI<3“TEQ'. 1") GOTO 605 

IF (KK3.EQ.NXT GO TO 777 ' 

WR IT E f6 79000) X S ( K 2T* A NS "" " 

ITTNPUNCHVE Q" • TTG’0T056 2 ' 

GOT 05 1 1 O' " " “ 

‘ 5'62 ' PUNCH9006 TXS (1<2') VANS 
GOT05TTO 
777 FANS2=ANS 


o> 

CD 
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GOT 051 1 0 

605_ ANS= (FANS3+ANS )_/2._ 

WRITE ('6 *9000 ) X S (~K 2 ) 7a N S _ ’ _ 

I F ( NPUNCH • EQ 1 ) GOT 0556 

G0T051 1 0 

"556 PUNCH9000 . XS (i<2 ) « ANS "" _ _ 

G0T051 1 0 """" ~~ " 

' " "778 IF (NCI .E'Q.'l ) GOTO 1228"" 

GOTO 1227 ^ 

1228 IF ( A3S (XS (<2 5-XR5 ) '• L E • .Ye-7')'g6t6"1 229 " 

I F ( ( NDX . EO • 0 • AND • XS (K2 ) . GT . XR5 )'. OR . ( NDX_» EQ - . 1 ."AND. XS’<K2 ) • LTV 

XX R 5 ) ) GOTO 1230 
GOTO 1 227 

1229 WR I TE (6 . 9000 )XS (K2 > • ANS 
IF (NPUNCH. EG. 1 VG0T0557 

1230 ANS=ANS/DX*DX5 ^ 

1227 WRITE(6.90C0 )XS CK2)".ANS 

\F (NPUNCH. EG. 1 )GOT0560 ' _ 

G0T05110 

557 PUNCH9000.XS (K2 ). ANS " "" " "" ' ' 

GOT 01230 

560 PUNCH9000 . XS (K2 ) , ANS 

GOT 051 1 0 " " " " 

779 IF ( KK3 •" EQ .1 ) GOT 06 1 1 1 

IF (NC2.EQ. 1 )G0T01 328 ' 

1327 WR I TE ( 6 '9000 )XS <K2 )", ANS 

I F (NPUNCH ."EG", f ) G0T0558 — _ 

GOT 051 10 ' .' ~ 

558 PUNCH9000 .XS < K2 ) . ANS 

gotosTTo 

6111 A N S = (FANS 2 + A NS ) /2 . 

GOTO 1327 

' 1 328 IF"( APS ( XS (K2 )"-X~R4 )". LEV. TE-77 GOT 01 329 

IF ( (NDX • EQ • 0 • AND • XS ( K 2 ) • GT » XR4 ') • OR . "( NDX . EQ . TTaND .~XsTl<~2T i 

XLT .XR4 iTGOfOl 330 * 

Gbt"bT327 

1329 WRITE (6. 9000 )XS"0<2) .ANS 

" IF (NPUNC'hVeQ . 1 ) G0T0559 

- 1 330 ANS="ANS7bX*bX4 
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782 WRITE(6*90C0 )XS(K2) *ANS 
j - ( NjD!j\icH'. EQ . 1 ) G0T056 1 

GO 70 5 1 10 

559 PUNCH9000 * XS ( l<2 ) * ANS 
GOT01330 

56 1 PUNCH9000 V XS (K2 ) *ANS~ 
5110' CONTINUE 
‘RETURN 

9000 FORMA 1 7~2 ET6Y8 f 
"‘END 


SUB R OUT I N E D I S C 0 T ( X A * Z A V TA dT* T AB Y if A S Z * N C 7NY « N Z VANS! “ 

Z T RE" dT MENS': IONS" 1 N~ Th rs“SUBROljT: I N E "AreTo nl y dummy dimensions. 

DIMENSION TABX (2 ) , TABY <2 ) * TABZ (2 ) *NPX (8 ) ,NPY (8 ) , YY (8 ) 

C DTmENSTON' TABX ( 2) ,'T ABY ( 2l *TABZ ( 2 ) * NPX ( 87 , N P yT6 “ ) ,7778“) 

calU“uns Tnc* i a' «7 dxY i'dzTT ms - : j 

IF "(NZ-17 " ' S ; 5 ♦ 1 0 " 

“5 CALL D'l SSER ( XA *TABX7l 77f*' N Y~, T DX * NNT 
NNN= I DX+ 1 

CALL' LAGRAN ( XA * TABX (NN ) * T ABY ( NN )* NNN , AnS V “ 

_ , - G0To 7Q - - - 

10~“ZARG = ZA ■’ """ " 

IP1X=IDX+1' 

TpTz^TdZ^T 
TF~TiXi 1~5 *~25 * 15 

1 5 'IF (ZARG-TABZTn'Z) V 25'*25i 20 

“20 ZARG"=TABZ ( NZ“ ) 

25” CALL “DT SSER - TZ‘ARG'» TABZ'( 'll *T ,~NZ“* Tl DZ7NPZ“> 

NX=NY/NZ 
NPZL=NPZ+ I DZ 

I F ' ( I MS )" 30*30*40 

30 CALL D I SSER (XA~* T ATbxTF f . T 7 N X",TD X 7 n P xTl ) ) 

D0~3fT" J J^NPZ * NPZL 

’npt? rr=uj-i“)TNx+NPX( tt ~ 

NPX ( j ) =NPxrn 
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-3 

CO 

35 "I = 1 + 1 ' 

Goto ' so 

4CTD0 ~4"5~ JUWPZVNPZL"”' ~ " " - - - - 

- rs =t j j- rr*Nx+i “ " 

CAUL 0 I SSFR ■' (XX, TABXrry , I S", NXlTDYiUPx ( 1 ) ) ' 

njpy'TT) =npx fTT ” ' " "" ' 

45 1 = 1 + 1 

-srarw 5 " 5 ' i = n rpyz - 

N[: - ocr 5 R px- m ~ - - - - 

NL0CY = NPVTT") 

' '55 -CALL" LA GRAN T XA fTABX' ('NUCC YVT'ABY ’( NLOC'Y' j i IP 1 XVYY U ) ) 
C AUC “ LAG RAN “ T"Z A R Gl TABZ {NP ZTY Y Y "( D Tl PTZI A NS j 
70 RETURN ~ 

’ END 


' " SUBROUTINE UNS T C »' I A W DX i I DZ , I MS ) ' 
~ I F _ (’ I C ) ■■ 515,^0 

5 I M S = 1 

""NC = -ic" "" 

goto'Ts 

ro""'i'Ms'=o'~ 

" NC=IC " 

15 IF "(NC-lOO) ' 20 * 25Y25”’" 

20 rA=o 

GOTO” 3‘0 

25 I A=1 ' - “ ' ■' ' - 

NC=NC-100 

30”TDX‘=NGZT"0‘ 

I DZ = NC- I DX* 1 0 

RETURN 

END 
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SUBROUTINE LAGRAN ( XA i X.< Y » N • ANS ) 
DIMENSION X(2)»Y(2) 

DIMENSION X ( 2 ) *Y(2) 

SUM = 0 *0 
DO 3 I = 1 « N 
PROD= Y ( I ) 

'DO" 2 J= 1 i N 
A = X ( I ) -X ( J ) 

IF (A ) 1 t2» 1 

1 B= (XA-X ( J) )/A 
PROD=PRQD*B 

2 CONTINUE 

3 SUM=SUM+PROD 
ANS=SUM 

'RETURN 

END 


SUBROUTINE DISSER ( XA i T AB » I »NXf I D « NPX ) 
DIMENSION TAB (2) 

D I MENS LON ' TAB (2 ) 

NPT = I D+ 1 

NPB = NPT/2 "" 

NPU=NPT-NPB 
IF ( NX-NPT ) ' 10*5 ( 10 

" 5 NPX= I ^ 

RETURN 

10 NLOW=I +NPB 

"J NUPP = I ^NX- (^PU+_1 ) 

DO 15'lT=NL0wTiNUPP ' 

_NLOC= 11 

IF (TAB(II)-XA) " 1 5 « 20 »20 
15 "CONTINUE” 

NPX=NUPP— NPB+ 1 

RETURN 

20 NL=NLOC-NPB 

NU=NL+ I D 


APPENDIX A 



-a 

^ _ 

"CO 25 J J = NL * NU 

NDIS=jj 

IF (TAB ( JJ )-TAB ( JJ+1 ) ) 25*30*25 

25“ _ C0NfiNuF 

________ - “ 

’ ’ RETURN ' ™ 1 " 

30 IF ( T AB ( N D I S )-XA )' 4_0^35r_35 __ 

"RETURN 
4 0 N P X = N D I S + 1 
-RETURN 
' END 
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FIXED INPUT DATA LISTING 



Gi 


FIXED INPUT DATA FOR SHARP CONES 

( 1 ) 

14 

NXL1 


6 

NTHN1 ( I ) , I = 1 , 1 4 


6 

o 

e 

7 

7 10 

1 1 

1 1 

12 

1 6 

1 6 

ML AM 1 ( I ) 

, I 

= 1,14 








35 

30 

48 

42 

42 60 

66 

66 

72 

96 

96 

MACH (1,1 

) * 

1 = 1,14 








.700 


.80 0 

.90 0 

.950 1 

. .000 

1.100 

1 

.200 

1 .500 

3.000 

n 

.500 4 

.000 

15. .coo 






XL 1 ( J ) , J 

= 1 

, 6 








.0 


. 2 

.4 

.6 

.3 

1 .0 




THN1 ( I , J ) , 

1 = 1,14 

J=1 « NTHN 1 ( I ) 






C. 0 


10.0 

15.0 

25.0 

33.0 

40.0 




c.o 


10.0 

20 .0 

3C .0 

40.0 





0.0 


10.0 

15.0 

20 .0 

25.0 

30.0 


33.0 

40.0 

0.0 


1 o.c 

15.0 

20 .0 

25.0 

33.0 


40.0 


0.0 


2.5 

10.0 

20.0 

25.0 

33.0 


40.0 


0.0 


2.5 

5.0 

7.5 

10.0 

20.0 


25.0 

30.0 

0.0 


2.5 

5.0 

7.5 

10.0 

15.0 


20.0 

25.0 

40.0 










o.c 


2.5 

5.0 

7.5 

10.0 

15.0 


20.0 

25,0 

40.0 










o.c 


2.5 

5 • C 

7.5 

10.0 

15.0 


20.0 

25.0 

35.0 


4 C . 0 








c.o 


2.5 

5.0 

7.5 

10.0 

12.5 


15.0 

17.5 

25.0 


27.5 

30.0 
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FIXED INPUT DATA FOR FRUSTUMS FOLLOWING CONE CYLINDERS 
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NTH N6 C I ) • I =1 f 8 

33333332 
NTHF6C I • J) t 1=1 *8 J = 1 » NTHN6 ( I ) 

4 4 4 

4 4 4 

4 4 4 

4 4 2 

2 4 2 

2 4 2 

2 4 2 

4 4 

NFA6 ( 1 i JtK ) i J=1 *3 K= 1 * NTHF6 ( 1 * J ) 

4 4 4 3 

4 4 4 4 

4 4 4 4 

NFA6 (2*J*K)*J=1*3 K = 1 ♦ NTHF6 ( 2 * J ) 

4 4 4 3 

4 4 4 4 

4 4 4 4 

NFA6 ( 3 * J * K ) * J = 1 t3 K= 1 ♦ N THF6 ( 3 ♦ J ) 

4 4 4 3 

4 4 4 4 

4 4 4 4 

NFA6<4* J*K ), J=1 ,3 K= 1 , NTHF6 ( 4 * J ) 

3 3 3 2 

4 4 4 4 

3 3 

NFA6(5% JtK ) * J=1 ,3 K= 1 * NTHF6 (5i J ) 

3 2 

4 4 4 4 

3 3 

NFA6(6*J*K ) *J = 1 *3 K = 1 « N THF6 ( 6 i J ) 

3 2 

4 4 4 4 
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NFA6 ( 7 , J 

iK),J 

= 1,3 

K= 1 , NTHF6 ( 7 , J ) 

3 

2 



4 

4 

4 

4 

3 

3 



NFA6C8, J 

, K ) , J 

= 1 ,2 

K=1 , NTHF6 ( 8 , J ) 

4 

4 

4 

4 

4 

4 

4 

4 

NLAM6 ( 1 , 

J ,K ) , 

J=1 , 

3 K = 1 , NTHF 6 ( 1 , J ) 

72 

72 

72 

54 

72 

72 

72 

72 

72 

72 

72 

72 

NLAM6<2, 

J,K ) , 

J= 1 , 

3 K = 1 ,NTHF6 (2,J ) 

72 

72 

72 

54 

72 

72 

72 

72 

72 

72 

72 

72 

NLAM6 ( 3 , 

J tK ) , 

J=1 i 

3 |< = 1 t NT HF 6 ( 3 , J ) 

72 

72 

72 

54 

72 

72 

72 

72 

72 

72 

72 

72 

NLAM6 ( 4 , 

JiK ) » 

J=1 , 

3 K=1 , NTHF 6 ( 4 « J ) 

54 

54 

54 

36 

72 

72 

72 

72 

54 

54 



NLAM6 ( 5, 

J , K ) i 

c 

n 

3 K = 1 , NTHF 6 ( 5 , J ) 

54 

36 



72 

72 

72 

72 

54 

54 



NLAM6 ( 6 , 

■ J t K ) , 

J=1 ♦ 

3 K = 1 , NTHF6 ( 6 , J ) 

54 

36 



72 

72 

72 

72 

54 

54 



KILAM6 ( 7 i 

iJiKli 

> J=1 ♦ 

3 K = 1 ♦ NTHF6 ( 7 ♦ J ) 

54 

36 



72 

72 

72 

72 

54 

54 



NLAM6 (8< 

» J , K ) i 

i J=1 i 

,2 K = 1 , NTHF 6 ( 8 , J ) 

72 

72 

72 

72 

72 

72 

72 

72 
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NT0TL6 (I»J).I=1»8 J=ltNTHN6(I) 


270 

288 

288 

270 

288 

288 

270 

288 

288 

1 98 

288 

108 

90 

288 

108 

90 

288 

108 

90 

288 

108 

288 

288 



NTOTFA ( I , J ) , I =1 ,8 J= 1 • NTHN6 ( I ) 


15 

1 6 

16 

15 

1 6 

16 

15 

16 

16 

1 1 

16 

6 

5 

1 6 

6 

5 

16 

6 

5 

16 

6 

16 

16 



MACH ( 6 . I ) . 1 = 1 <8 

.80 .90 1.20 1.50 2.18 


XD6 ( J ) * J = 

1*18 



o 

• 

o 

.2 

• 4 

.6 .8 

3.00 

3.50 

4.00 

4.50 5.00 

THN6 ( I i J ) 

*1=1*8 

J=1 » NTHN6 ( I ) 

1 5.0 

22.5 

30.0 


o 

• 

in 

22.5 

30.0 


15.0 

22.5 

30.0 


1 5 • 0 

22.5 

30.0 


15.0 

22.5 

30.0 


15.0 

22.5 

30.0 


15.0 

22.5 

30.0 


15.00 

22.50 



THF 6 ( 1 * J* 

K ) * J = 1 

*3 K= 1 

* NTHF6 ( 1 * J ) 

5.0 

10.0 

15.0 

20.0 

5.0 

10.0 

15.0 

20.0 

5.0 

10.0 

15.0 

20.0 

THF6 (2. J. 

K ) * J=1 

*3 K= 1 

* NTHF6 ( 2 * J ) 

5.0 

10.0 

15.0 

20.0 

5.0 

10.0 

15.0 

20.0 

5.0 

10.0 

15.0 

20.0 


<o 

-j 


2.81 

4.04 

15.00 



1 .0 
5.50 

1 .25 
6.00 

1.50 

6.50 

2.00 

2.50 
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THF6(3, J 

• K ) , J= 1 

iH 

II 

y 

n 

*■ 

i NT HF 6 (3 

. J ) 

o 

• 

in 

10.0 

15.0 

20.0 


in 

• 

o 

10.0 

15.0 

20.0 


5.0 

10.0 

15.0 

20.0 


T HF 6(4. J 

• K ) t J = 1 

.3 K= 1 

t NTHF 6 ( 4 

^ J > 

5.0 

10.0 

15.0 

20.0 


5.0 

10.0 

15.0 

20.0 


o 

• 

in 

20.0 




THF6 ( 5 * J 

« K ) . J = 1 

.3 K = 1 

.NTHF 6 ( 5 

. J ) 

5.0 

20.0 




5.0 

10.0 

15.0 

20.0 


01 

• 

o 

20.0 




THF 6 ( 6 * J 

. K ) . J= 1 

.3 K = 1 

.NTHF6 ( 6 

. J ) 

o 

• 

in 

20.0 




01 

• 

o 

10.0 

15.0 

20.0 


o 

• 

in 

20.0 




THF6(7, J 

* K ) * J = 1 

.3 K= 1 

.NTHF6 (7 

. J ) 

o 

• 

in 

20.0 




o 

• 

in 

10.0 

15.0 

20.0 


o 

• 

in 

20.0 




THF6 ( 8 * J 

» K ) . J = 1 

.2 K = 1 

.NTHF6 (8 

« J ) 

5 • 00 

10.00 

15.00 

20.00 


5.00 

10.00 

15.00 

20.00 



FA6 ( 1 « J. 

K ) 

. J=1 .3 

K = 1 .1 

NTOTFA ( 1 t J ) 








0.0 


1 .0 

2.0 

4.0 

0.0 

1 .0 

2.0 

4.0 

o 

• 

o 

1 

• 0 

2.0 


4.0 

1 .0 

2.0 

4.0 







o 

. 

o 


1 .0 

o 

. 

CM 

4.0 

0.0 

1 .0 

2.0 

4.0 

o 

. 

o 

1 

• 0 

2.0 


4.0 

0.0 

1 .0 

2.0 

4.0 






0.0 


1 .0 

2.0 

4.0 

0.0 

1 .0 

o 

. 

OJ 

4.0 

0*0 

1 

• 0 

2.0 


4.0 

O 

. 

o 

1 .0 

O 

. 

CVJ 

4.0 






FA6 (2. J. 

K ) 

. J=1 .3 

K= 1 * I 

NTOTFA ( 2 i J ) 








0.0 


1 .0 

2.0 

4.0 

o 

. 

o 

1 .0 

o 

. 

CM 

4.0 

o 

. 

o 

1 

• 0 

o 

. 

C\J 


4.0 

1 .0 

2.0 

4.0 







o 

• 

o 


1 .0 

2.0 

4.0 

o.o 

1 .0 

O 

. 

CM 

4.0 

o 

. 

o 

1 

• 0 

o 

. 

0J 


4.0 

0.0 

1 .0 

2.0 

4.0 






0.0 


1 .0 

2.0 

4.0 

0.0 

1 .0 

O 

. 

CM 

4.0 

o 

• 

o 

1 

• 0 

2.0 


4.0 

o 

. 

o 

1 .0 

O 

. 

CVI 

4.0 






FA6 (3. J. 

K) 

. J = 1 . 3 

K = 1 . 

NTOTFA (3, J ) 








0.0 


1 .0 

2.0 

4.0 

o 

. 

o 

1 .0 

2.0 

4.0 

o 

. 

o 

1 

• 0 

o 

. 

CM 


4.0 

1 .0 

2.0 

4.0 
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0.0 

1 .0 

2.0 

4.0 0.0 

2.0 

4.0 

0.0 

1.0 2.0 

0.0 

1 .0 

2.0 

4.0 0.0 

2.0 

4.0 

0.0 

1.0 2.0 

FA 6 ( 4 . J . K ) 

» J = N 3 

K =1 * 

NTOTFA (4 .J ) 

o 

» 

o 

1.0 

4.0 

0.0 1.0 

4.0 




0.0 

1.0 

2.0 

4.0 0.0 

2.0 

4.0 

0.0 

1.0 2.0 

0.0 

1 .0 

4.0 

0.0 1.0 

FA 6 ( 5 . J . K ) 

* J = 1 *3 

K = 1 ♦ 

NTOTFA ( 5 . J ) 

o 

• 

o 

1 .0 

4.0 

1.0 4.0 

0.0 

1 .0 

2.0 

4.0 0.0 

2.0 

4 . C 

0.0 

1.0 2.0 

C.O 

1.0 

4.0 

0.0 1.0 

F A 6 ( 6 « J * K ) 

i J =1 #3 

K >1 • 

NTOTFA ( 6 * J ) 

0.0 

1.0 

4.0 

1.0 4.0 

0.0 

1 .0 

2.0 

4.0 0.0 

2.0 

4.0 

0.0 

1.0 2.0 

0.0 

1 .0 

4.0 

0.0 1.0 

FA 6 ( 7 . J . K ) 

• J = 1 .3 

K = 1 .1 

NTOTFA ( 7 . J ) 

0.0 

1.0 

4.0 

o 

. 

o 

. 

M 

0*0 

1 .0 

2.0 

4.0 O.o 

2.0 

4.0 

0.0 

1.0 2.0 

0.0 

1.0 

4.0 

0.0 1.0 

FA 6 ( 8 * J « K ) 

t J =1 .2 

K = 1 . 1 

NTOTFA ( 8 . J ) 

0.0 

1 .0 

2.0 

4.0 0.0 

2.0 

4.0 

0.0 

1.0 2.0 

0.0 

1.0 

2.0 

-P> 

• 

O 

o 

• 

o 

2.0 

4.0 

0.0 

1.0 2.0 


K 1 ( 1 , J , K ) • J =1 ,3 K = 1 . NTOTFA ( 1 , J ) 

0 0 1 1 0 0 1 

0 

0 0 1 1 0 0 0 

0 0 

0 0 0 0 0 0 0 

0 0 

K 1 ( 2 , J . K ) , J = 1 ,3 K = 1 . NTOTFA ( 2 . J ) 

0 0 0 0 0 0 

0 


CD 

CD 


0 


1.0 

4.0 

2.0 

4.0 

0.0 

1 .0 

1.0 

4.0 

2.0 

4.0 

O 

• 

O 

1 .0 

4.0 

0.0 

1 .0 

4.0 

1 .0 

1 .0 
4.0 
4.0 

O 

. 

OJ 

4.0 

0.0 

1.0 

1.0 

4.0 

4.0 

o 

. 

0J 

4.0 

O 

• 

o 

1 .0 

1.0 

4.0 

4.0 

2.0 

4.0 

o 

• 

o 

1.0 

1.0 

4.0 

4.0 

o 

. 

OJ 

4.0 

o 

. 

o 

1 .0 

1 .0 
4.0 

o 

. 

CVJ 

4.0 

o 

. 

o 

1 .0 

1 .0 

o 

. 

OJ 

4.0 

o 

• 

o 

1.0 


1 0 0 0 0 0 1 
1 0 0 0 0 0 0 
0 0 0 0 0 0 0 

0 0 0 0 0 0 0 
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100 


0 0 0 0 0 0 0 

0 0 

0 0 0 0 0 0 0 

0 c 

K1 (3, JiK ) , J=1 ,3 K=1 ,NT0TFA(3» J) 

0 0 0 0 0 0 0 

0 

0 0 0 0 0 0 0 

0 0 

0 0 0 0 0 0 0 

0 0 

K1 (4i JtK ) ij=l « 3 K = 1 » NTOTF A ( 4 » J ) 

0 110 110 

0 1110 11 

1 1 

0 110 11 

K1 (5*J*K)«J=1 i3 K = 1 ,NT0TFA(5» J) 

0 1111 

C 1 1 1 0 1 1 

1 1 

0 110 11 

K1 (6* JiK ) * J=1 i3 K = 1 ,NTOTFA (6i J ) 

0 1111 

0 1110 11 

1 1 

0 110 11 

K1 (7* J»K ) « J=1 ,3 K=1 ,NT0TFA<7» J) 

0 1111 

0 1110 11 

1 1 

0 110 11 

K1(8»J*K),J=1,2 K=1 ,NT0TFA(8, J) 


0 

1 

1 

1 

0 1 

1 

1 

1 





0 

1 

1 

1 

0 1 

1 

1 

1 





L AM6 ( 1 » 1 

*K ) 

*K = 1 

.270 



1 • 10 


.23 

• 50 

1 .00 

.97 

.53 


• 63 

.79 

• 38 

• 50 

• 69 


.43 

.46 

• 52 

• 64 


ooooooo 

0 0 0 0 0 0 0 

ooooooo 

ooooooo 

ooooooo 

1111 
10 1110 1 

10 1110 1 

10 1110 1 

10 1110 1 


1 

0 

1 

1 

1 

0 

1 

1 

0 

1 

1 

1 

0 

1 


CO 

CM 

• 

• 53 

• 74 

• 87 

.34 

h* 

If) 

. 

.74 

CVJ 

• 

• 47 

.53 

.44 

• 47 

• 53 

• 62 

.48 
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• 51 

.54 

.58 

.60 

.57 

• 56 

.63 

.60 

.58 

. 56 

.52 

.58 

.56 

.59 

• 52 

.56 

.59 

.52 

.58 

.56 

.52 

.58 

1 .33 

.38 

1 .03 

1 .24 

.63 

1.12 

1 .22 

1 .23 

1 .12 

1.10 

1.18 

.84 

1.12 

1.13 

.94 

1 .22 

1 .07 

1.10 

1 .24 

1 .23 

1 • 1 0 

1.27 

1 .20 

1.10 

1 . 2.2 

1 .20 

1 .23 

1.10 

1 .20 

1 .23 

1.10 

1 .22 

1 .20 

1.10 

1 .22 

1 .20 

1 .23 

1.41 

1.13 

1 .46 

1 .41 

.77 

1.18 

1 .41 

.96 

1 .22 

1.27 

1 .41 

i — > 

• 

ro 

m 

1 .23 

1 .41 

1.18 

1 .28 

1 .4 1 

1 .34 

1 .39 

1 .36 

1 .41 

1 .31 

1 .38 

1 .41 

1 .36 

1 .33 

1 .4 1 

1 .36 

1 .33 

1 .38 

1 .41 

1 .33 

1 .38 

1 .41 

1 .36 

1 .33 

1 .20 

1 .60 

.96 

1 .25 

1 .52 

1 .43 

1 .26 

1 .30 

1 .42 

1 .34 

1 *46 

1 .50 

1 .46 

1 .50 

1 .51 

1 .49 

1.51 

1 .50 

1 .49 

1.51 

1 .51 

1 .50 

1 .49 

1 .51 

1 .50 

L AM 6 ( 1 * 2 « 

< ) .K=l 

9288 



1 .20 

.06 

.95 

1.10 

1 .06 

• 53 

.64 

.86 

.4 1 

.48 

.73 

.52 

.56 

.52 

.68 

. 62 

.56 

.59 

• 60 

• 62 

.52 

• 60 

.63 

.54 

.52 

• 63 

.54 

• 52 

.60 

.63 

.52 

. 60 

.63 

.54 

.52 

. 63 

. 54 

1 .44 

.20 

1 .40 

1 .29 

.72 

1.13 

1.13 

1 .24 

1 .08 

1 .00 

1.18 

.96 

1 .08 

1.13 

.99 

1.10 

1 .03 

1.10 

1.19 

1.17 

1.10 

1 .02 

1 .20 


.59 

.56 

. 62 

.60 

.60 

.63 

.57 

• 58 

.56 

.59 

.58 

.56 

• 59 

.52 

.58 

.59 

.52 

• 58 

.56 

.59 

1 .50 

1 .29 

.51 

1.10 

1 .33 

.73 

1.13 

1.13 

1 .20 

o 

CO 

. 

1.08 

1.14 

.88 

1.17 

1.06 

1.12 

1 .26 

1.14 

1 .08 

1.24 

1 .23 

1 • 1 0 

1 .27 

1.18 

1.23 

1.22 

1 .20 

1 .23 

1.10 

1.22 

1 .23 

1.10 

1 .22 

1 .20 

1 .23 

.47 

1 .07 

1 .60 

1 .41 

• 63 

1 .36 

1 .41 

• 88 

1 .21 

1.30 

1 .03 

1 • 22 

1 .25 

1 .41 

1.12 

1 .23 

1 .41 

1 .26 

1 .38 

1.31 

1 .36 

1 .34 

1 .38 

1 .41 

1.36 

1 .38 

1 .41 

1 .36 

1 .33 

1.38 

1 .36 

1 .33 

1 .38 

1 .41 

1.36 

1 .38 

• 62 

1.10 

1 .70 

.81 

1 .08 

1 • 27 

1 .46 

1.19 

1 .28 

1 .36 

1 .41 

1 .49 

1 .42 

1.42 

1.50 

1 *50 

1 .49 

1 .51 

1 .50 

1 .50 

1 .49 

1 .51 

1 .50 

1.49 

1 .49 

1 .51 

1 .50 

1 .49 

1 .51 

.22 

. 66 

• 78 

.94 

.33 

.57 

• 78 

.47 

.52 

.53 

.56 

.61 

.53 

• 64 

.60 

.60 

.54 

• 60 

• 61 

• 60 

.60 

• 63 

.54 

.52 

.60 

.54 

• 52 

.60 

• 63 

.54 

.60 

. 63 

.54 

.52 

• 60 

1 .63 

1 .37 

.54 

1 .23 

1.30 

. 84 

1*10 

1 .04 

1 .21 

.92 

1 .00 

1.14 

.98 

1 .09 

1.0 1 

1 .00 

1.15 

1.10 

1 .08 

1.00 

1 .20 

1.10 

1 .02 

. 

ro 

o 

1.20 
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1 • 10 

1 .02 

1 .20 

1 .20 

1.10 

OJ 

o 

• 

1 .20 

1 .20 

1.10 

1 .02 

1 .20 

1 .20 

1.10 

1 .02 

1 .20 

1 .20 

1.10 

1 .02 

1 .20 

1 .20 

1.10 

1 .02 

1 .20 

1 .20 

1 .67 

.31 

1 .56 

1 .80 

1 .62 

.77 

1 .43 

1 .50 

1 .57 

.98 

1 .33 

1 .36 

1 .54 

1.10 

1 .28 

1 .27 

1.51 

1.18 

1 .26 

1 .22 

1 .48 

1 .23 

1 .27 

1 .20 

1 .45 

1 .25 

1 .28 

1 .20 

1 .43 

1 .27 

1 .32 

1 .22 

1 .39 

1 .27 

1 .37 

1.31 

1 .38 

1 .27 

1 .40 

1 .37 

1 .40 

1 .27 

1 .40 

1 .40 

1 .40 

1 .27 

1 .40 

1 .40 

1 .40 

1 .27 

1 .40 

1 .40 

1 .40 

1 .27 

1 .40 

1.40 

1 .40 

1 .27 

1 .40 

1 .40 

1 .40 

1 .27 

1 .40 

1 .40 

1 .40 

1.27 

1 .40 

1 .40 

1 .40 

1 .27 

1 .40 

1 .40 

2.00 

.51 

1.72 

1.94 

1 .97 

1 .03 

1 .59 

1 .67 

1 .94 

1 .27 

1 .51 

1 .53 

1 .91 

1.38 

1 .47 

1.46 

1 .89 

1 .46 

1 .47 

1 .42 

1 .87 

1 .50 

1 .48 

1 .40 

1 .84 

1 .53 

1 .50 

1 • 42 

1 .83 

1 .54 

1 .53 

1 .43 

1 .81 

1 .56 

1 .58 

1 .51 

1 .80 

1 .57 

1 .60 

1 .57 

1 .80 

1 .58 

1 .56 

1.60 

1 .80 

1 .58 

1 .56 

1 .60 

1 .80 

1 .58 

1 .56 

1 .60 

1 .80 

1.58 

1 .56 

1 .60 

1 .80 

1 .58 

1 .56 

1 .60 

1 • 80 

1 .58 

1 .56 

1 .60 

1 .80 

1 .58 

1 .56 

1 .60 

1 .80 

1 .58 

1 .56 

1 .60 



LAM 6 ( 1*3 

* K ) * K = 1 

*288 








1 .52 

0.00 

.91 

1 .02 

1 .22 

*29 

.71 

.83 

1.16 

• 38 

.61 

.73 

.94 

• 44 

• 56 

♦ 67 

• 87 

.47 

• 54 

.64 

.83 

• 50 

.54 

• 63 

• 80 

• 53 

• 54 

• 63 

• 80 

• 57 

.56 

.63 

.83 

• 60 

• 59 

.64 

• 94 

.58 

• 60 

.63 

1 .06 

.56 

• 57 

• 57 

1 • 06 

• 56 

• 57 

• 57 

1 • 06 

.56 

.57 

.57 

1 • 06 

• 56 

• 57 

.57 

1 .06 

• 56 

• 57 

• 57 

1 .06 

• 56 

• 57 

• 57 

1 .06 

.56 

• 57 

.57 

1 .06 

• 56 

.57 

.57 

1 *86 

• 00 

1 .34 

1 .48 

1 • 67 

• 81 

1.19 

1 .29 

1 .54 

.99 

1.10 

1 • 1 8 
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1 .92 

2.02 

1 .93 

1 .83 

1 .92 

2.10 

1 .95 

1 .79 

1 .92 

2 . 10 

1 .95 

1.79 

1 .92 

2. 10 

1 .95 

1 .79 

1 .92 

2 . 1 0 

1 .95 

1 .79 



L AM 6 ( 6.3 

,K ) ,K = 1 

,108 








.52 

.84 

.20 

• 52 

• 50 

.13 

• 52 

.36 

• 12 

.52 

.34 

. 1 2 

.52 

.35 

• 13 

.52 

.38 

.14 

.52 

.40 

• 16 

• 52 

.42 

.18 

• 52 

.46 

.24 

.52 

.48 

.32 

.52 

.49 

.40 

.52 

• 50 

.47 

.52 

.50 

.51 

.52 

.48 

.55 

.52 

.47 

• 57 

.52 

.42 

.55 

.52 

.39 

.45 

.52 

.39 

.45 

1 .92 

1 .42 

1 .51 

1 .92 

1 .44 

1 .40 

1 .92 

1 .48 

1 .33 

1 .92 

1 .50 

1 .30 

1 .92 

1 .53 

1 .30 

1.92 

1 .58 

1.31 

1 .92 

1 .61 

1 .36 

1 .92 

1 .65 

1 .41 

1 .92 

1.72 

1 .55 

1 .92 

1 .80 

1 .68 

1 .92 

1 .88 

1 .78 

1 .92 

1 .93 

1.82 

1 .92 

1 .99 

1 .87 

1 .92 

2.03 

1 .83 

1 .92 

2.09 

1 .79 

1 *92 

2.09 

1 .79 

1 .92 

2.09 

1 .79 

1 .92 

2.09 

1 .79 



L AM 6 ( 7.1 

iK ) ,K = 1 

,90 








.52 

.74 

.30 

.52 

.72 

.30 

.52 

.70 

• 30 

.52 

. 68 

.30 

.52 

.65 

. 30 

.52 

.62 

• 30 

• 52 

.60 

.30 

.52 

.58 

.30 

.52 

.53 

.32 

• 52 

• 49 

.36 

.52 

.47 

.38 

.52 

.43 

.40 

. 52 

.42 

• 41 

.52 

• 45 

.42 

.52 

.48 

.43 

.52 

. 54 

.43 

• 52 

• 60 

• 45 

.52 

.60 

.45 

2.05 

1 .80 

1 .90 

1 .60 

1.80 

1.49 

1 .76 

1 .42 

1 .72 

1 .39 

1 .73 

1 .38 

1 .76 

1 .41 

1 .81 

1 .48 

1 .95 

1 .64 

2.08 

1 .80 

2 . 19 

1 .92 

2.25 

1 .99 

2.28 

2.01 

2.30 

1 .98 

2.30 

1 .90 

2.30 

1 .90 

2.30 

1.90 

2.30 

1 .90 
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L AM 6 ( 7,2 

* K > * K = 1 

,288 



.52 

.72 

• 40 

.30 

.52 

.42 

.30 

.52 

.65 

.43 

.52 

.61 

.44 

.30 

• 52 

• 45 

• 30 

.52 

.54 

.44 

.52 

.50 

.41 

• 38 

• 52 

.4 1 

.41 

.52 

• 46 

• 42 

.52 

.43 

.50 

.43 

.52 

• 53 

.45 

1 .02 

1 .22 

1 .08 

1 .02 

1 .03 

.82 

• 58 

1 .02 

.70 

.50 

1 .02 

.88 

• 66 

1 .02 

.81 

.62 

.50 

1 .02 

.73 

.66 

1 .02 

.79 

• 83 

1 .02 

.87 

1 .00 

.91 

1 .02 

1 .02 

.98 

1 .02 

1 .08 

.97 

1 .02 

1 .09 

.95 

.94 

1 .52 

1 .33 

.98 

1 .52 

1 .33 

1 .20 

1 .50 

1 .20 

1 .05 

.92 

1 .52 

1 .02 

.95 

1 .52 

1.19 

1 .04 

1 .52 

1.31 

1 .30 

1.19 

1 .52 

1 .54 

1 .39 

1 .52 

1 .58 

1 .59 

1 .52 

1 .83 

1 .45 

1 .46 

1 .52 

1 .36 

1 .43 

1 .52 

1 .90 

1 .36 

1 .98 

1 .73 

1 .77 

1 .60 

1 .98 

1 .53 

1 .42 

1 .98 

1 .58 

1.49 

1 .98 

1 .63 

1 .50 

1 .41 

1 .98 

1 .72 

1 .64 

1 .96 

2.10 

1 .97 

1 .98 

2.24 

2.33 

1 .99 

1 .98 

2.39 

1 .98 

1 .98 

2.27 

2.32 

1 .98 

2.27 

2.32 

1 .90 

1 .98 

LAM 6 ( 7,3 

• K ) ,K = 1 

, 1 08 



• 52 

.51 

.30 

• 52 

. 56 

• 48 

• 30 

.52 

• 44 

• 30 

.30 

• 52 

.38 

.30 

• 52 

.52 

.37 

.38 

.52 

. 38 

.43 

.42 

.52 

.48 

.43 

.45 

.52 

.57 

.45 

1 .98 

1 .98 

1 .50 

1 .49 

1 .98 

1 .43 

1 .4 1 

1 .38 

1 .98 

1 .47 

1.41 

1 .64 

1 .98 

2.02 

1 .80 

1 .98 


70 

.41 

• 30 

• 52 

• 68 

30 

.52 

.62 

• 44 

• 30 

60 

. 45 

• 30 

.52 

.58 

32 

• 52 

.51 

.43 

.36 

48 

• 40 

• 40 

.52 

.47 

42 

• 52 

.44 

• 48 

• 43 

43 

.53 

.45 

• 52 

• 43 

73 

1 .02 

1 .12 

.93 

• 64 

97 

.75 

.52 

1 .02 

.92 

49 

1 .02 

• 84 

.63 

.48 

78 

• 65 

.57 

1 .02 

.78 

76 

1 .02 

.82 

• 92 

• 84 

92 

1 .03 

.98 

1 .02 

1.00 

96 

1 .02 

1 .09 

• 95 

.94 

60 

1 .60 

1 .03 

1 .52 

1 .45 

94 

1 .52 

1 .25 

1.10 

.93 

18 

1 .02 

.93 

1 .52 

1.18 

98 

1 .52 

1 .23 

1.15 

1 .08 

38 

1 .43 

1 .30 

1 .52 

1 .48 

46 

1 .52 

1 .71 

1 .57 

1 .48 

90 

1 .36 

1 .43 

1 .52 

1.90 

43 

1 .98 

1 .88 

2.03 

1.80 

64 

1 .62 

1 .49 

1 .98 

1 .60 

39 

1 .98 

1 .58 

1 .48 

1.38 

72 

1 .54 

1 .48 

1 .98 

1 .93 

80 

1 .98 

2. 19 

2.08 

1 .92 

27 

2.40 

2.01 

1 .98 

2.27 

90 

1 .98 

2.27 

2.32 

1 .90 

27 

2.32 

1 .90 



> 30 

.52 

.51 

• 30 

CVJ 

if ) 

. 

.52 

• 42 

.30 

• 52 

.40 

37 

• 32 

• 52 

• 36 

.36 

40 

.52 

.40 

.41 

.52 

52 

. 52 

.43 

• 52 

.57 

81 

1 .80 

1 .98 

1 . 60 

1.60 

42 

1 .96 

1 .40 

1 .39 

1 .98 

98 

1 .55 

1 .48 

1 .98 

1 .78 

20 

1 .92 

1 .98 

2.31 

1 .99 
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1 .98 

2.40 

2.01 

1 .98 

00 

• 

CM 

1 .98 

1 .98 

2.48 

1 .90 


1.98 

2.48 

1 .90 

1 .98 

2.48 

1.91 

1 .98 

2.48 

1 .91 




L AM 6 ( 8*1 

• K ) « K =1 

,288 









.543 

.561 

.561 

.561 

.543 

.557 

• 557 

• 557 

.543 


.555 

.555 

.555 

.543 

.553 

.553 

.553 

.543 

• 552 

• 552 


.552 

.543 

.551 

.551 

.55 1 

.543 

.549 

.549 

• 549 

.543 


• 548 

• 548 

.548 

.543 

. 546 

.546 

.546 

.543 

• 545 

.545 


• 545 

• 543 

.545 

.545 

.545 

.543 

.545 

. 545 

.545 

• 543 


• 544 

• 544 

. 544 

.543 

.543 

.543 

.543 

.543 

.543 

.543 


• 543 

• 543 

.543 

.543 

.543 

.543 

.543 

.543 

.543 

.543 


.543 

.543 

.543 

1 .070 

1 .333 

1 .333 

1 .333 

1 .070 

1 .298 

1 .298 

1 

.298 

1 .070 

1 .275 

1 .275 

1 .275 

1 .070 

1 .257 

1 .257 

1 .257 

1 .070 

1 

.243 

1 .243 

1 .243 

1 .070 

1.070 

1 .070 

1 .070 

1.070 

1 .070 

1 .222 

1 

• 070 

1 .070 

1 .070 

1 .208 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 • 070 

1 

• 070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 

• 070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 

• 070 

1 .070 

1 .070 

1 . 070 

1 .070 

1 .070 

1 .070 

1 .070 

1.070 

1 .070 

1 

• 070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .564 

2.408 

2.408 

2.408 

1 .564 

2 

.294 

2«294 

2.294 

1 .564 

2.205 

2.205 

2.205 

1 .564 

1 .564 

2.144 

2 

• 144 

1 .564 

1 .564 

2.097 

2.097 

1 .564 

1 .564 

1 .564 

2.051 

1 .564 

1 

• 564 

1 .564 

2.010 

1 .564 

1 .564 

1 .564 

1 .966 

1 .564 

1 .564 

1 .564 

1 

• 564 

1 .564 

1 .564 

1.564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 

.564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 

• 564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 . 564 

1 .564 

1 .564 

1 

• 564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

2.011 

3.753 

3.753 

3 

• 753 

2.01 1 

3.459 

3.459 

3.459 

2.01 1 

2.011 

3 % 269 

3.269 

2.011 

2 

• 01 1 

3 . .137 

3.137 

2.011 

2.011 

2.011 

3.014 

2.011 

2.01 1 

2.011 

2 

• 917 

2.011 

2.011 

2.01 1 

2.401 

2.011 

2.011 

2.01 1 

2.011 

2.01 1 

2 

.01 1 

2.011 

2.011 

2.011 

2.011 

2.011 

2.01 1 

2.011 

2.01 1 

2.01 1 

2 

• 01 1 

2.01 1 

2.01 1 

2.01 1 

2.01 1 

2.01 1 

2.011 

2.011 

2.01 1 

2.01 1 

2 

.01 1 

2.011 

2.011 

2.01 1 

2.01 1 

2.011 

2.011 

2.011 

2.01 1 

2.01 1 

2 

• 01 1 

2.011 

2.011 

2.01 1 

2.011 

2.01 1 

2.011 

2.011 

2.01 1 
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,K ) *K = 1 
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.543 

.561 

.561 

.561 

.543 

.557 

.557 

.557 

.543 


.555 

.555 

.555 

.543 

.553 

.553 

.553 

.543 

• 552 

• 552 


• 552 

.543 

.551 

.551 

.551 

• 543 

.549 

.549 

.549 

.543 


• 548 

.548 

.548 

.543 

.546 

.546 

• 546 

.543 

.545 

.545 


• 545 

.543 

.545 

.545 

• 545 

.543 

.545 

.545 

.545 

.543 


• 544 

.544 

. 544 

.543 

.543 

.543 

.543 

• 543 

.543 

.543 


.543 

.543 

.543 

.543 

.543 

.543 

.543 

.543 

.543 

.543 


.543 
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.543 

1 .070 

1 .333 

1 .333 

1 .333 

1 .070 

1 .298 

1 .298 

1 .298 

1 .070 

1 .275 

1 .275 

1 .275 

1 .070 

1 .257 

1 .257 

1 .257 

1 .070 

1 .243 

1 .243 

1 .243 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 *070 

1 .222 

1 .070 

1 .070 

1 .070 

1 .208 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 #070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1.070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .070 

1 .564 

2.408 

2.408 

2.408 

1 .564 

2.294 

2.294 

2.294 

1 .564 

2.205 

2.205 

2.205 

1 .564 

1 .564 

2.144 

2.144 

1 .564 

1 .564 

2.097 

2.097 

1 .564 

1 .564 

1 .564 

2.051 

1 .564 

1 .564 

1 .564 

2.010 

1 .564 

1 .564 

1 .564 

1 .966 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 .564 

1 # 56.4 

1 .564 

1 .564 

1 .564 

1 .564 

2.011 

3.753 

3.753 

3.753 

2.0 11 

3.459 

3.459 

3.459 

2.01 1 

2.011 

3.269 

3.269 

2.01 1 

2.011 

3.137 

3.137 

2.011 

2.011 

2.011 

3.014 

2.011 

2.011 

2.01 1 

2.917 

2.011 

2.01 1 

2.01 1 

2.401 

2.0 1 1 

2.011 

2.011 

2.0 11 

2.01 1 

2.01 1 

2.011 

2,01 1 

2.011 

2.01 1 

2.01 1 

2.01 1 

2*011 

2.011 

2.01 1 

2.011 

2.011 

2.01 1 

2.011 

2.011 

2.01 1 

2.011 

2.011 

2.011 

2.01 1 

2.011 

2.011 

2.011 

2.011 

2.011 

2.01 1 

2.0 1 1 

2.011 

2.01 1 

2.01 1 

2.011 

2.011 

2.01 1 

2.01 1 

2.01 1 

2.01 1 

2.011 

2.011 

2.011 
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FIXED INPUT DATA FOR CYLINDERS FOLLOWING FRUSTUMS 

NMACH ( 7 ) 

8 

NXD7 

24 

NDDR7 ( I ) * I = 1 ♦ 8 

1 1 1 3 3 3 3 2 

NTHF7C 1 « J) « 1 = 1 * 8 J=1»NDDR7(I) 

4 ‘ 

4 

4 

2 2 2 

3 3 2 

3 3 2 

3 3 2 

2 2 

NLAM7< I . J) . 1 = 1 .8 J=1 » NDDR7 ( I ) 

96 

96 

96 

48 48 48 

72 72 48 

72 72 48 

72 72 48 

48 48 

MACH ( 7 . I ) * 1 = 1 <8 


o 

00 

• 

1 .00 

1 .20 

1 .50 

2.50 

3.50 

5.00 

1-5.00 



XD7 ( 1 ) * I 

= 1.24 








0.0 

• 1 

.2 

.3 

.4 

.5 

• 6 

.7 

• 8 

.9 

1 .00 

1 .25 

1 .50 

1.75 

2.00 

2.50 

3.00 

3.50 

4.00 

4.50 


5.00 5.50 6.00 6.50 

DDR7( I . J) . 1 = 1 .8 J=1 i NDDR7 ( I ) 


.70 



.70 



.70 



.25 

.50 

.75 

in 

OJ 

• 

• 50 

.75 

in 

OJ 

• 

0 

in 

• 

.75 

in 

OJ 

. 

0 

in 

. 

.75 

.25 

0 

in 

. 
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THF7 < 1 » J»K > « J=1 » NDDR7 ( 1 ) K= 1 » NTHF 7 ( 1 ♦ J ) 

5.0 10.0 15.0 30.0 

THF7 (2 . J iK ) . J=1 .NDDR7 (2 ) K= 1 « NTHF7 ( 2 » J ) 

5.0 10.0 15.0 30.0 

THF7 ( 3 « J i K ) » J=1 » NDDR7 ( 3 ) K=1 i NTHF7 ( 3 * J ) 

5.0 10.0 15.0 30.0 

THF7 ( 4 « J i K ) . J=1 . NDDR7 ( 1+ ) K= 1 i NTHF7 (4i J) 

5.0 10.0 

5.0 10.0 

5.0 10.0 

THF7 ( 5 » J » K ) , J = 1 . NDDR7 ( 5) K= 1 , NTHF 7 (5i J ) 

5.0 10.0 20.0 

5.0 10.0 20.0 

5.0 20.0 

THF7(6«JiK ) i J=1 . NDDR7 (6) K= 1 » NTHF7 ( 6 ♦ J ) 

5.0 10.0 20.0 

5.0 10.0 20.0 

5.0 20.0 


THF7 ( 7 »_J (K ) «_J= 1 NDDR7 ( 7 ) K=_l •NTHF 7 < 7 ♦ J ) 

5.0 10.0 20. 0 

5.0 10.0 20.0 

5.0 20.0 

T HF7 (8.J.K ) »J=1 ,NDDR7(8) K= 1 « NTHF7 ( 8 « J ) 

5.00 10 . 00 . 

5.00 10.00 


LAM7 ( 1 . 1 »K ) »K=1 .96 


- .370 -1.140 

-1.140 

o 

0 
10 

• 

f—t 

1 

-.070 

.120 

%570 

-1 .400 

%04Q 

.160 

.580 -1 .300 

.075 

. 145 

.400 

-.220 

.087 

%\ 30 

.205 

.600 

. 088 .110 

• 1 60 

.685 

.079 

. 100 

. 128 

.590 

*062 

• 080 

.100 #400 

. 055 

• 070 

.080 

.310 

%050 

%060 

« 072 

• 245 

•048 .058 

.070 

.200 

.050 

.055 

.0 65 

• 1 55 

• 060 

• 060 

.060 .130 

.070 

.050 

.050 

• 120 

.070 

.085 

.090 

.147 

.060 .088 

.092 

. 1 30 

.059 

.080 

.080 

• 120 

%060 

• 070 

.070 .128 

.060 

.075 

.072 

.120 

• 060 

• 080 

.080 

.110 

.060 .087 

.080 

.095 

.060 

.087 

.080 

.095 

%060 

• 087 

.080 .095 

• 060 

.087 

.080 

.095 





LAM 7 (2*1 «K ) »K=1 

.96 








.347 1.200 

.510 

.1 30 

.345 

.630 

%360 

0 
r- 

. 

1 

%340 

.460 

.320 -.140 

. 330 

.620 

.530 

.340 

%300 

%620 

.740 

.818 

.030 .540 

.750 

.855 

-.060 

.350 

.540 

.840 

-.118 

• 140 
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.240 

.720 

0 
co 

. 

1 

-.020 

.060 

.280 

-.215 

-•160 

-.060 

• 040 

-.210 

-.248 

-. 170 

-.090 

-.110 

-.280 

-.320 

-.240 

-.030 

-.255 

-.333 

-.30 0 

.020 

-. 120 

-.210 

-. 190 

.050 

• 080 

-.1 13 

-.070 

.070 

.120 

-.080 

.033 

.080 

.120 

-.025 

• 080 

• 085 

.115 

.040 

.115 

.085 

.100 

.070 

.132 

• 080 

.092 

• 068 

.140 

.068 

.090 

••060 

.125 

• 055 

. 090 

.060 

.100 

• 055 

• 090 

.060 

.100 

.055 

.090 

.060 

. 100 





L AM 7 ( 3 » 1 

iK )» K=l 

% 96 








.580 

.720 

.450 

.260 

.470 

.500 

.200 

-.440 

• 410 

• 420 

.150 

-.430 

.380 

.460 

.270 

• 450 

.340 

.458 

.350 

.520 

.295 

.430 

.390 

.390 

.235 

.390 

.430 

• 452 

• 180 

.350 

.410 

.450 

.130 

.300 

.355 

.440 

.090 

.240 

.290 

.390 

.058 

• 175 

.240 

.300 

-.0 10 

.040 

• 130 

.145 

-.020 

• 000 

.060 

.070 

.015 

-.0 10 

.035 

.040 

-.015 

-.020 

• 020 

.020 

.010 

- .01 5 

.010 

-.010 

• 060 

.020 

.050 

• 005 

• 1 10 

• 070 

. 1 00 

.045 

. 135 

.110 

• i 45 " 

• 100 

.1 10 

• 130 

• 159 

.130 

.020 

.050 

• 100 

.060 

• 000 

• 050 

.040 

• 050 

%000 

• 050 

.040 

.050 

.000 

.050 

• 040 

• 050 





L AM 7 ( 4 . 1 

«K ) tK = l 

% 48 








.374 

.723 

♦ 355 

. 652 

• 338 

• 590 ' 

' . 320 

• 534 

• 304 

.485 

.286 

.440 

.272 

• 406 

• 258 

• 372 

' .243 

• 340 

• 226 

.317 

.214 

• 290 

.181 

.242 

• 145 

• 208 

• 128 

• 173 

• 1 10 

*146 

.087 

.118 

• 046 

.059 

%030 

• 033 

.018 

.018 

• 010 

• 010 

.007 

.007 

.005 

.005 

• 003 

• 003 

• 003 

%003 
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» K ) iK=l 

.48 








.249 

.710 

.251 

.620 

• 253 

• 548 

.250 

• 484 

• 248 

.430 

.244 

.387 

• 240 

.352 

• 237 

.318 

.230 

.290 

.222 

.263 

<215 

.240 

. 1 95 

. 197 

. 1 70 

• 170 

• 145 

• 145 

• 130 

.130 

.115 

.115 

.084 

.084 

.070 

.070 

.058 

• 058 

• 045 

.045 

.036 

.036 

.025 

-025 

.0 15 

.015 

.0 1 5 

.015 



LAM 7 ( 4*3 

• K ) ♦ K = 1 

448 








.220 

.220 

.208 

.208 

. 1 94 

• 1 94 

.185 

.185 

• 174 

.174 

• 1 65 

. 1 65 

.157 

. 157 

• 1 47 

.147 

. 140 

.140 

• 133 

*1 33 

.128 

• 128 

.115 

.115 

.100 

.100 

.0 92 

.092 

.086 

*086 

• 077 

.077 

.060 

.060 

.052 

.052 

.045 

.045 

.040 

• 040 

.037 

.037 

.035 

.035 

.030 

.030 

.030 

.030 



L AM 7 ( 5 * 1 

♦ K ) * K =1 

.72 








.277 

.530 

.805 

• 274 

.503 

.738 

.272 

.474 

.678 

.270 

.450 

.626 

.267 

.415 

.576 

.264 

• 403 

• 534 

• 261 

.382 

.500 

.257 

.364 

.470 

.255 

.375 

.440 

• 252 

.329 

.415 
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.246 

.313 

.390 

.235 

. 1 86 

.345 

.224 

.260 

• 310 

• 210 

.240 

.276 

• 194 

.218 

• 250 

• 159 

• 181 

.201 

.122 

• 149 

• 160 

.088 

. 120 

. 122 

.060 

.092 

.092 

.040 

• 070 

• 070 

.026 

.054 

. 054 

• 015 

.041 

.045 

• 008 

.034 

.034 

• 008 

.034 

.034 
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.233 

.486 

.561 

• 234 

.4 66 

.538 

• 234 

.446 

♦ 512 

.233 

.426 

.490 

.231 

.408 

• 464 

.229 

.385 

.438 

.226 

.366 

.4 1 5 

.225 

• 350 

.393 

.220 

.330 

.372 

.217 

.313 

• 350 

.210 

.284 

.330 

.200 

.256 

.282 

.184 

.220 

• 242 

• 164 

.190 

.205 

• 1 49 

. 162 

• 1 74 

.117 

. 127 

.128 

• 094 

• 100 

• 1 00 

.080 

.080 

.080 

.065 

• 065 

.065 

• 050 

• 050 

• 050 

.038 

.038 

• 038 

.028 

.028 

.028 

.015 

.015 

.015 

• 015 

.015 

.015 
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• 1 80 

.180 

• 1 73 

• 173 

• 1 66 

.166 

• 1 60 

• 160 

• 1 54 

• 154 

.138 

. 1 38 

• 144 

. 144 

.138 

• 138 

.133 

• 133 

.127 

• 127 

• 124 

.124 

.112 

• 112 

• 104 

• 1 04 

.096 

.096 

.080 

• 080 

.079 

.079 

.069 

.069 

.062 

.062 

• 060 

.060 

.055 

.055 

.052 

• 052 

.052 

.052 

• 050 

• 050 

• 050 

.050 
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• 278 

• 540 

.820 

.276 

• 510 

• 750 

.274 

.480 

.688 

• 272 

• 452 

.632 

• 270 

.428 

• 586 

• 266 

• 405 

.542 

.264 

• 385 

• 508 

• 260 # 

.366 

.476 

.258 

.346 

.446 

.254 

.33 1 

.422 

.250 

. 315 * 

.396 

.239 

.286 

• 350 

.228 

• 260 

• 314 

• 213 

.238 

.280 

• 1 97 

.220 

.250 

.161 

. 185 

• 2 06 

.120 

• 150 

• 1 43 

.090 

• 120 

. 128 

• 065 

.095 

• 0 95 

• 045 

.072 

• 072 

.030 

.055 

• 055 

.019 

• 045 

.045 

.010 

• 032 

• 032 

• 010 

• 032 

.032 
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.254 

.412 

.490 

.250 

.400 

.472 

.246 

.389 

.455 

.239 

• 372 

.440 

.232 

• 360 

• 424 

.228 

• 348 

• 4 08 

.220 

• 335 

.390 

.216 

.323 

.377 

• 2 1 0 

.310 

.360 

.205 

.299 

• 345 

.199 

.286 

.330 

• 186 

.257 

.291 

. 172 

.230 

• 260 

• 160 

.203 

.230 

• 150 

. 180 

• 1 96 

• 128 

.139 

. 1 48 

. 1 08 

• 1 13 

• 1 15 

.090 

.090 

.090 

• 072 

.072 

.072 

.058 

.058 

• 058 

.042 

.042 

.042 

.028 

• 028 

.028 

.015 

• 015 

.015 

• 015 

• 015 

.015 
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. 1 35 

. 1 33 

. 133 

• 131 

.131 

.129 

.129 

.127 

• 127 
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.125 

.125 

.121 

.121 

.119 

.119 

.118 

.1 18 

.115 

• 1 15 

• 112 

.112 

.108 

. 108 

. 1 04 

.104 

.100 

.100 

.092 

• 092 

• 064 

.084 

• 075 .. 

.075 

.067 

.067 

.060 

•060 

• 059 
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• 054 
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.055 
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.054 

.054 
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.• 340 . . 
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.500 
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.346 

.466 
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.332 

.434 

• 31 1 
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.310 

.376 

.290 

.298 
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.280 

.289 
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.273 

• 277 

.312 

.262 

• 265 

.294 

.256 

.260 

.280 

.246 
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.268 

.231 

.233 

.240 
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.220 

.220 

.200 

.216 
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.209 
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.•187 
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.170 

. 170 
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.092 
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.070 

• 070 

• 050 

• 070 
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• 220 

.261 

.290 

.220 

.260 

.287 

.220 

.257 

.280 

.220 

• 252 

.272 

.220 

.243 

.267 

.220 

.239 

.260 

.219 

.234 

.253 

.216 

.231 

.250 

• 216 

.216 

.245 

.213 

.220 

.239 

• 210 

.208 

.232 

.204 

.204 

.220 

.195 
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.208 
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• 1 86 
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.063 
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.045 
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.028 

.028 
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• 146 

• 146 

.143 
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.133 

.133 

• 131 
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• 1 30 

• 130 

• 1 28 

.128 

. 126 

.126 

.123 
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• 1 19 

.119 

.113 

.113 

. 1 08 

.108 

• 1 02 

.102 
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.072 

.072 
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.060 

. 058 
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.054 

.054 

.056 

.056 

.054 

.054 
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.116 

.116 

.116 

.116 

.116 
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.116 
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.116 

.116 

.116 

.116 

.116 

.116 

.116 

.116 

.116 
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.116 

.116 

.116 

.116 

.116 

. 1 16 

. 1 16 
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. 1 16 
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.116 
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• 1 16 
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• 1 16 

.116 

.116 
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.116 

.116 

.116 

.116 

.116 

.116 

.116 

.116 

.116 

.116 

.116 

. 1 16 
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FIXED INPUT DATA FOR BOATTAILS FOLLOWING LONG CYLINDERS 


NMACH (8 ) 

a 










o 

NXD 8 










12 










NTHB 8 ( I ) 

* 1 = 1*8 









5 

4 5 

4 

5 5 

4 

5 





NLAM 8 ( I ) 

* 1 = 1*8 









60 

48 60 

48 

60 60 

48 

60 





MACH ( 8 * I 

) * 1 = 1 *8 









1 • 50 

1.91 

2.00 

2.50 

3.00 

5.00 

6.00 

15.00 



XD 8 < I ) * I 

= 1*12 









• 00 

• 10 

*20 

.40 

• 60 

.80 

1 .00 

1 .25 

1.50 

1.75 

2.00 

2.50 









THB 8 ( IiJ ) iI = l »8 . 

J = 1 » NTHB 8 ( I ) 







0.0 

4.0 

8.0 

12.0 

16.0 






0.0 

5.6 

7.0 

9.3 







o.o 

4.0 

8.0 

12.0 

16.0 






0.0 

4.0 

8.0 

12.0 







0.0 

2.0 

4.0 

8.0 

12.0 






0.0 

2.0 

4.0 

8.0 

12.0 






0.0 

6.0 

12.0 

18.0 







0.0 

4.0 

8.0 

12.0 

16.0 






LAM 8 ( 1 ♦ J ) , J =1 *60 









.000 

.000 

-.030 

-.035 

-.060 

.000 

-•020 

-.060 

-.085 

-•040 

.000 

-.048 

-. 105 

-. 150 

-•032 

• 000 

-.120 

-.240 

-.320 

-.120 

.000 

-.220 

-.412 

— . 660 — 

1 .200 

• 000 

-•261 

-.590 

- 1.100 

- 1 .480 

.000 

-.250 

-.790 

- 1 .330 

• 000 

• 000 

-.170 

-.860 

- 1 .400 

-.240 

.000 

-.100 

-.820 

- 1.380 

-• 1 50 

• 000 

-.060 

-.760 

-1 .300 

--•090 

.000 

-.040 

-.710 

- 1.190 

-.060 

• 000 

-.040 

-.710 

- 1.190 

-.060 

L AM 8 ( 2 t J ) » J = 1 .48 









.000 

-.050 

-. 1 55 

-. 157 

• 000 

-.080 

-.168 

-.170 

.000 

-.1 12 

-.182 

-.189 

.000 

-.174 

-.220 

-.230 

.000 

-.235 

-.263 

-.283 

• 000 

-.295 

-.318 

-.355 

• 000 

-.346 

-.350 

-.393 

.000 

-.378 

-.354 

-.410 

.000 

-.380 

-.354 

-.420 

.000 

-.369 

-.352 

-.417 

.000 

-.358 

-.348 

-.396 

.000 

-.358 

-.348 

-.396 



L AM 8 ( 3 ♦ J ) iJ = 1 .60 









.000 

-.030 

-.030 

-.090 

• 000 

.000 

-.048 

-.062 

-.132 

-•048 

.000 

-.064 

-.095 

-. 180 

-• 105 

• 000 

-.100 

-.160 

-.280 

-.240 

.000 

-.129 

-.232 

-.380 

-.390 

.000 

-.155 

-.312 

-.500 

-.470 
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• 000 

1 

• 

1 —* 

-J 

Ul 

1 

.408 

-.630 

-•170 

• 000 

-.160 

-.450 

-•650 

-•090 

.000 

-.110 - 

.450 

-.620 

-•050 

.000 

-.060 

-.410 

-.570 

-•020 

• 000 

-.030 - 

.400 

-.530 

-.0 10 

• coo 

-.030 

-•400 

-.530 

-•010 

LAMB ( 4 * J 

) * J =1 *48 









.000 

-.040 

. 075 

-.110 

.000 

-.050 

-.100 

-.140 

• 000 

-•065 

-.120 

-. 170 

.000 

-.092 

-.1 59 

-.228 

.000 

-.116 

-.1 88 

-.280 

• 000 

-.130 - 

.215 

-.325 

• 000 

-•130 

-.240 

-•368 

• 000 

-• 120 

-.240 

-.390 

.000 

-#100 

-.210 

-.390 

.000 

-.060 

-.160 

-•340 

.000 

-.030 - 

. 130 

-.290 

• 000 

-.030 

-.130 

-.290 



L AM 8 ( 5 * J 

) * J = 1 * 60 









• 000 

.000 

.000 

.000 

.000 

.000 

-.010 

-.020 

- 030 

- 050 

.000 

-.020 - 

. 040 

-.079 

-.1 20 

.000 

-.038 

-.050 

- 165 

- 280 

• 000 

-.040 

. 1 30 

-.280 

-.470 

.000 

-.048 

-• 1 80 

- 418 

- 720 

• 000 

-.058 - 

.225 

-.557 

- 1 .020 

• 000 

-.023 

-•300 

-.760 

-1 • 480 

• 000 

-.065 ~ 

.390 

-.940 

- 2.020 

.000 

-.065 

-.450 

- 1.140 

- 2.150 

.000 

-.065 - 

.500 

- 1 .300 

- 2 . 1 50 

.000 

-.065 

-.500 

- 1 .300 

- 2.150 

LAM 8 ( 6 . J 

) * J =1 *60 









.000 

.000 

.000 

.000 

• 000 

• 000 

.000 

-.010 

-.010 

-•010 

.000 

-.020 - 

.030 

-.032 

-.042 

.000 

-.050 

-.062 

-.080 

-.1 15 

.000 

-.075 - 

.090 

-. 1 33 

-.198 

• 000 

-.0 98 

-.1 08 

-.200 

-•325 

.000 

-.105 - 

.135 

-.282 

-.550 

• 000 

-. 120 

-•180 

-.410 

-.81 0 

.000 

-.150 - 

.222 

- . 545 

-.950 

• 000 

-.130 

-.268 

-.730 

-.930 

.000 

-.100 - 

.319 

-•920 

-.875 

• 000 

-. 100 

— • 31*9 

-.920 

-•875 

LAMS ( 7 . J 

) * J =1 *48 









• COG 

-.034 - 

. 075 

-.250 

.000 

-.125 

-.4 13 

-.612 

• 000 

-.108 

- *346 

-.529 

.000 

-.080 

-.245 

-.383 

.000 

-.062 

-.180 

-.303 

.000 

-.053 - 

.131 

-.240 

.000 

-.049 

-.088 

-•182 

• 000 

-.049 

-.050 

-.100 

.000 

-.049 

-.050 

-.050 

.000 

-.049 

-•050 

-•050 

.000 

-.049 - 

.050 

-.050 

• 000 

-.049 

-.050 

-•050 



LAM 8 ( 8 * J 

) * J = 1 *60 









.000 

• 003 

.000 

• 000 

• 000 

• 000 

.003 

• 000 

• 000 

• 00.0 

• 000 

.003 

.000 

.000 

• 000 

.000 

.003 

• 000 

• 000 

• 000 

.000 

.003 

.000 

.000 

.000 

.000 

.003 

• 000 

• 000 

• 000 

.000 

.003 

.000 

.000 

.000 

.000 

.003 

.000 

• 000 

• 000 

.000 

.003 

.000 

.000 

.000 

.000 

• 003 

.000 

• 000 

• 000 

.000 

.003 

.000 

.000 

.000 

.000 

.003 

.000 

• 000 

• 000 
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